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ABSTRACT 


Electrical  networks  consisting  of  lumped  linear  and  memoryless  non- 
linear  elements  and  an  arlitrary  number  of  lossless  transmission  lines 
are  considered  It  is  shown  that  a  large  class  of  such  networks  may  be 
described  by  a  system  of  functional-differential  equations  having  the  form 

x(t)  >e  f(xt) 

where  the  state  cf  the  system  at  time  t  >  0  is  represented  by  x^,  a  point 

in  the  space  Cu(  (-<*>,  0],  E11  )  cf  bounded  continuous  functions  mapping  the 
a 

interval  (-oo,cJ  into  En  ,  with  the  compact  open  topology,  and  the  function 

f  mapping  C,.(  {•».  0],  En  )  into  En  :is  continuous  and  locally  Lipschitz ian . 
n 

A  Lyapunov  f'uncticnal  is  presented  and  use'  t,  obtain  several  theorems  con¬ 
cerning  the  stability  and  instability  cl'  the  equilibrium  solution,  x  =  0, 
cf  the  network  Several  examples  the  theory  are  presented. 
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Chapter  I 


INTRODUCTION 

In  this  dissertation  several  theorems  are  presented  which  state 
sufficient  conditions  to  ensure  that  an  equilibrium  state  of  a  nonlinear 
distributed  network  is  stable,  asymptotically  stable,  completely  stable, 
or  unstable.  We  use  the  name  "nonlinear  distributed  network"  to  refer 
to  an  electrical  network  consisting  of  lumped  linear  and  memoryless 
nonlinear  elements  and  an  arbitrary  number  of  lossless  transmission 
lines.  The  technical  literature  abounds  with  such  theory  for  lumped 
nonlinear  networks  [1,2];  while,  for  linear  networks,  both  lumped  and 
distributed,  the  problem  is  simply  that  of  locating  the  roots  of  the 
network's  characteristic  equation'1'.  To  the  author's  best  knowledge, 
relatively  little  has  been  written  concerning  the  more  general  case 
of  nonlinear  distributed  networks. 

There  appear  to  be  at  least  two  ways  of  obtaining  stability 
criteria  for  nonlinear  distributed  networks:  One  way  is  to  write  the 
partial  differential  equations  which  govern  the  distributed  elements 
of  the  network  and  then  consider  as  boundary  conditions  or  constraints , 
the  algebraic  and  ordinary  differential  equations  which  arise  by 


Much  has  been  written  on  methods  for  determining  the  location  of 
the  zeros  of  exponential  sums,  that  is,  functions  of  the  complex  variable 

n  C^.z 

z  of  the  form:  =  2  A^(z)e  ,  where  the  A^  are  polynomials  in  z, 

k=0 

and  the  are  constants.  The  characteristic  equations  for  linear 
networks  containing  lo..  .ess  transmission  lines  are  of  this  form.  The 
reader  is  referred  to  references  [3.>k,5,6],  and  especially  Chapters  12 
and  13  of  reference  [7],  where  further  references  are  to  be  found.  For 
jertain  kinds  of  linear  distributed  networks,  several  authors  have 
obtained  sufficient  conditions  to  ensure  that  all  I'oots  of  the  charac¬ 
teristic  equation  have  negative  real  parts.  See,  for  example,  [8]. 
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applying  Kirchoff'c  laws  to  the  lumped  portion  of  the  network.  Then, 
applying  the  Lyapunov  theory  for  dynamical  systems  [9]  (of  which  the 
boundary  value  problem  is  a  particular  kind)  stability  criteria  may 
be  obtained.  This  approach  has  been  considered  by  Brayton  and 
Miranker  [10].  On  the  other  hand,  one  may  treat  each  distributed 
element  in  the  network  as  a  two -port  and  obtain  mathematical  expressions 
which  show  the  manner  in  which  the  electrical  variables  at  the  ports  are 
related.  These  relations,  along  with  the  current-voltage  relations  of 
the  lumped  elements,  may  then  be  introduced  into  the  Kirchoff's 
voltage  and  current  law  equations  to  obtain  a  system  of  functional- 
differential  equations  which  describes  the  behavior  of  the  network 
(functional-differential  equations  are  obtained  because  the  expressions 
which  relate  the  electrical  variables  at  the  ports  of  the  distributed 
elements  are  functional  equations).  The  Lyapunov  theory  for  functional- 
differential  equations  may  then  be  applied  to  this  system  to  obtain 
stability  criteria.  This  second  approach  is  the  one  that  we  shall 
consider. 

Functional-differential  equations,  and  the  application  of  Lyapunov' 
second  method  for  determining  the  stability  of  solutions  of  these  equa¬ 
tions,  have  been  treated  by  several  authors  in  the  mathematical  liter¬ 
ature  [11 ,12,1.3] .  Recently,  J.K.  Hale  has  published  several  theorems 
[12,13]  which  we  have  found  to  be  particularly  suitable  for  the  kind 
of  functional- differentiae  equations  which  describe  a  large  class  of 
nonlinear  distributed  networks.  We  make  use  of  three  of  these  theorems 
in  our  work  and  will,  state  them  in  the  next  chapter. 

Let  us  now  consider  a  simple  example  which  will,  serve  to  demon 
strata  the  purpose  ana  scope  of  this  work.  We  start  with  the  lumped 
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network  of  figure  1.1,  which  *z  shall,  at  first,  assume  t.c  be  linear. 

We  also  assume  H,C  >  0.  The  resistor  r  any  be  negative.  Obviously, 
v  *  0  is  an  equilibrium  state  for  this  network.  If  we  wish  to  determine 

'{ 

Figure  1.1.  A  simple  lumped  linear  network. 

whether  or  not  this  is  a  stable  equilibrium,  the  procedure  in  very 
simple:  we  compute  the  value  of  conductance  for  the  parallel  combina¬ 
tion  of  resistors;  and.  If  it  is  positive,  the  equilibrium  is  stable. 
Letting  g  3  l/r  and  G  »  l/R  ve  have  g||G  ■  g  +  G  and  hence  our  stability 
criterion  is:  If  g  >  -G,  the  equilibrium  state  is  stable. 

Let  us  now  *dd  a  lossless  transmission  line  to  our  network,  as  in 
,  Ct,  and  i  denote,  respectively,  the  inductance 
per  unit  length,  the  capacitance  per  unit  length,  and  the  length  of 
the  transmission  line.  With  this  modification  we  may  find  that  a 

a 

*(s) 

b 

Figure  l.c.  A  simple  distributed  linear  network. 


Figure  1.2  where 


3 


stable  equilibria  state  has  become  unstable.  For  example,  let 
r  -  -1.78,  R  -  1.42,  C  »  1/9,  Lt  «  81/37,  Ct  -  37/81,  *  «  s/6.  We  see 
that  r  <  -R,  which  implies  that  g  >  -0,  and  therefore  the  equilibrium 
state  v  -  0  tor  the  lumped  network  of  Figure  1.1  is  stable.  Clearly 
v  a  0  is  also  an  equilibrium  state  for  the  distributed  network.  We 
shall  now  show  that  this  equilibrium  state  is  unstable.  Since  the 
network  is  linear,  what  must  be  shown  is  that  there  exists  at  least 
one  root  of  the  network's  characteristic  equation  which  has  a  positive 
real  part.  It  is  a  well-known  fact  [14  pp.90-94,  15  pp. 262 -264]  that  the 
zeros  of  the  admittance  Y(s)  seen  at  the  port  a-b  in  Figure  1.2  are 
roots  of  this  network' 3  characteristic  equation.  In  fact,  due  to  the 
particularly  simple  nature  of  this  network,  these  are  the  only  roots 
of  the  characteristic  equation.  The  admittance  Y(s)  is  given  by  the 
formula 


X(aJ 


r 


sC 


(R  +  Z0)eTS  -  (R  -  Z^e”*8 
_(R  +  Z0,eT3  ♦  (R  -  7.o)e~*a 


i 


where  ZQ  =  »  2-19  ond  t  =  a  I  -  V  i  =  g  •  Now,  Y(s) 

has  zeros  at  approximately  s 


hjis 


(t  ±  -ip 


-ih  ±  if) 


(rf  1  if)  -(fi  +  if) 

3.eie  U  6  +0.77e  12  6 


=  -  0.561  +  0.056  +  jo. Ill  +  0.457[  1.104  +  JO. 241] 

-  (-0.561  +  0.056  +  0.505}  +  J(0.U1  -  o.ll') 

0 

Thus,  for  uie  distributed  netwarx,  the  equilibrium  state  v  n  0  is  unstable- 
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It  seed  not  always  happen  that  the  addition  of  lossless  trans¬ 
mission  lines  to  a  lumped  network  with  a  stable  equilibrium  state  causes 
the  resulting  distributed  network  to  have  an  unstable  equilibrium,  in 
particular,  if  c  ■  1  in  our  example,  both  the  lumped  and  distributed 
networks  have  stable  equilibria  v  »  0.  It  is  obvious  that  changing 
the  capacitor’s  value  has  no  effect  on  the  stability  of  the  lumped  net¬ 
work.  For  the  distributed  network,  however,  all  of  the  zeros  of  Y(b) 
now  have  negative  real  parts.  This  fact  is  proved  in  Appendix  A. 

It  is  interesting  to  note  that  if  the  lumped  network  of  Figure  1.1 
has  an  unstable  equilibrium  then  the  distributed  network  of  Figure  1.2 
also  has  an  unstable  equilibrium.  That  is,  the  addition  of  a  lossless 
transmission  line  to  our  lumped  network,  when  it  is  unstable,  cannot 
make  it  stable.  This  fact  Is  proved  in  Appendix  A. 

Hie  stability  criterion  for  the  linear  lumped  network  of  Figure  1.1 
may  be  expressed  graphically  as  In  Figure  1.3a.  A  straight  line  is 
drawn  in  the  i-v  plane  (vher  i  ind  v  denote,  respectively,  the  current 


Figure  l.J.  Stability  criteria  for  lumped  and  distributed  networks. 
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through  and  the  voltage  across  the  resistor  r)  whose  location  is  deter¬ 
mined  by  the  parameter  R  (G  --  l/p).  This  line  and  the  line  v  =-•  0  divide 
the  plane  into  two  regions.  If  the  i-v  curve  of  the  resistor  r  (which 
in  our  case  is  the  straight  line  i  =  gv)  lies  in  the  region  which  is 
labeled  "stable",  then  the  network  of  Figure  1.1  has  a  stable  equilibrium 
state  v  =  0.  If  the  i-v  curve  of  the  resistor  r  lies  in  the  region 
labeled  "unstable",  then  the  network  of  Figure  1.1  has  an  unstable  equi¬ 
librium  state. 

Some  of  the  stability  criteria  wh?ch  we  shall  derive  will  be  similar 

to  this  simple  criterion.  When  our  theory  is  applied  to  the  distributed 

network  of  Figure  1.2,  two  lines  are  determined  in  the  i-v  plane.  The 

position  of  these  lines  depends  only  on  the  parameters  R  and  Z  .  These 

o 

lines,  together  with  the  line  v  =  0,  divide  the  plane  into  three  regions 
as  shown  in  Figure  1.3b.  We  allow  the  resistor  r  to  have,  in  fact, 
almost  any  reasonable  (nonlinear)  i-v  curve.  Our  results  are  (in  part): 

1)  If  the  i-v  curve  for  the  resistor  r  lies  in  the  region  labeled 

"stable"  in  Figure  1.3b,  then  the  equilibrium  state  v  -  0  for  the  distrib- 

1 

uted  network  of  Figure  1.2  is  completely  stable  . 

2)  If  the  i-v  curve  for  the  resistor  r  lies  in  the  region  labeled 
"unstable"  in  Figure  1.3b,  then  the  equilibrium  state  v  -  0  for  the  distrib¬ 
uted  network  of  Figure  1.2  is  unstable. 

3)  If  the  i-v  curve  of  the  resistor  r  lies  in  the  remaining  region  of 
Figure  1.3b,  then  it  is  uncertain  whether  or  not  the  equilibrium  state 
v  =  0  for  the  distributed  network  of  Figure  1.2  is  stable, 

^Complete  stability,  sometimes  referred  to  as  asymptotic  stability 
in  the  large,  is  asymptotic  stability  where  the  region  of  asymptotic 
stability  is  comprised  of  all.  the  points  from  which  a  motion,  or  trajectory, 
may  originate  [16  p.  8,  1.7  pp. 56-66]. 
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Several  important  features  of  our  results  should  be  emphasised  at 
this  point.  The  following  remarks  refer  particularly  to  our  example} 
they  apply ,  however,  for  the  most  part,  to  all  networks  which  sure  members 
of  the  rather  large  class  for  which  the  theory  is  applicable.  First  of 
ail,  since  our  system  is  nonlinear,  it  is  significant  that  the  theory  can 
guarantee  complete  stability.  If  the  system  were  linear  it  would  be 
satisfactory  to  guarantee  only  asymptotic  stability  since  linearity 

4 

would  then  imply  complete  stability.  Since  all  systems,  however,  are 
to  some  extent  nonlinear,  it  is  really  complete  stability  (or  at  least 
asymptotic  stability  with  some  knowledge  of  the  extent  of  asymptotic 
stability)  that  is  needed  for  practical  applications  [17  P-57)-  Next, 
we  should  not  be  too  surprised  to  find  that  there  exists  a  region  in 
the  i-v  plane  for  which  the  stability  question  ia  left  unresolved.  The 
criteria  are  based  upon  a  knowledge  of  only  the  values  of  the  resistor 
R  and  the  characteristic  impedance  ZQ.  In  the  numerical  example  it  was 
found  that,  with  all  other  parameters  fixed,  a  network  could  be  made 
stable  or  unstable  by  adjusting  the  value  of  the  capacitor  C.  It  is, 
therefore,  rather  surprising  that  regions  can  indeed  be  found  where 
complete  stability  and  instability  can  be  assured  regardless  of  the 
value  of  so  many  parameters,  e.g.,  C,  a,  /.  On  the  other  hand,  it  is 
clear  that  if  all  of  the  parameter  values  were  considered  the  stability 
criteria  that  one  might  develop  would  not-  be  nearly  so  easy  to  apply  as 
our  results.  In  many  instances  it  might  be  significant  that  our  criteria 
do  not  depend  upon  the  length  of  the  transmission  line-  For  instance, 
if  one  were  designing  circuits  which  would  be  interconnected  by  trans¬ 
mission  lines,  it  might  be  important  to  know  that  the  resulting  network 
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would  be  stable  regardless  o*  the  actual  length  of  these  lines.  We  should 
also  point  out  at  thi3  time  that  the  stability  and  instability  regions 
obtained  by  our  methods  are  not  necessarily  the  best  regions  which  might 
be  found.  For  example,  it  is  obvious  from  physical  considerations  that 
the  region  labeled  "stable"  in  Figure  1.3b  should  always  contain  the 
first  and  third  quadrants.  For  certain  values  of  R  and  ZQ,  however,  the 
line  i  *  k^{R,Z  )v  has  positive  slope.  Tftis  unfortunate  circumstance 
does  not  occur  in  what  is  hoped  will  be  the  more  usual  situations;  that 
is,  when  the  transmission  line  is  terminated  in  a  resistor  R  whose  value 
is  somewhat  close  to  the  characteristic  impedance  of  the  line.  If,  in 
fact,  0-5  <  R/Zq  <  2.0  it  turns  out  that  the  stable  region  will  include 
the  first  and  third  quadrants.  Another  characteristic  of  our  results  is 
that  the  regions  labeled  "stable"  and  "unstable"  in  Figure  1.3b  are  always 
contained  in  the  corresponding  regions  in  Figure  1.3a.  This,  of  course, 
is  to  be  expected  since,  as  was  pointed  out  earlier,  the  theory  does  not 
take  into  account  the  length  of  the  transmission  line.  If  transmission 
lines  of  infinitesimal  length  were  present  in  a  network  of  lumped  elements, 
it  is  obvious  from  physical  considerations  that  the  behavior  of  the  net¬ 
work  should  approximate  that  of  the  corresponding  lumped  network.  If, 
in  our  example,  P/ZQ  1.  the  regions  labeled  "stable"  and  "unstable" 
in  Figure  1.3b  approach  the  corresponding  regions  in  Figure  1.3a.  This 
is  a  satisfying  result  since,  when  ZQ  ~  R,  the  distributed  network  is 
equivalent,  for  most  purposes,  to  the  lumped  network.  Thus,  our  results 
limbered  1)  and  c,'  above  apply  also  to  the  lumped  network  If  the  references 
tc  Figures  1  }d  and  1  j  .are  changed  to  read  Figures  l-3a  and  1.1  respec- 
1 ively  Since  v.;  have  shewn  (in  Appendix  A)  that  when  the  networks  oi 
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Figures  1.1  and  1.2  are  lineur,  the  distributed  network  is  unstable 
whenever  the  lumped  network  In  unstable,  one  might  expect  the  s nan  to 
apply  an  well  when  r  is  a  nonlinear  resistor.  Our  results  bear  this 
out  whenever  R  <  ZQ,  however,  we  do  not  guarantee  this  property  (although 
it  might  ati.ll  be  true}  when  R  > 
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Chapter  II 

LYAPUNOV  STABILITY  THEORY  FOR 
FIMCT 10 HAL  -VI FFEREHTIAL  EQUATIONS 

In  thia  chapter  we  give  precise  definitions  of  the  terms  which  will 
he  used  in  our  application  of  the  stability  theory  of  functional -differen¬ 
tial  equations.  Host  of  the  terms  used  are  standard  ones  in  the  mathe¬ 
matical  literature;  see,  for  example,  [11,18,19]-  We  also  state  three 
theorems  Of  J.K.  Hale  [12,13],  upon  which  much  of  our  work  is  based.  For 
proof  of  these  theorems  the  reader  is  referred  to  Hale's  paper, 
reference  [12]. 

1.  Basic  Definitions  and  Notation 

The  real  n-dimensional  Euclidean  jpace  is  denoted  by  E° ,  and  |lx||  ■ 
(xj[  +  ...  +  denotes  the  norm  of  an  element  x  in  E°.  The  elements 

x  in  E°  are  taken  to  be  column  vectors  and  denotes  the  corresponding 
row  vector.  Similarly,  if  M  is  an  m  ..  n  matrix  mapping  En  into  e”  then 
m'“  denotes  the  transpose  of  M;  also,  ||M||  denotes  the  norm  of  M,  defined 
by  JlMp  =  sup(jlHx|l:  x  <_  En,  |lxjl  «  1  ).  With  the  above  definition  for 
jjx  |; ,  it  can  be  shown  [20  pp.  59-60]  that  |!M|!  ^  vx  where  X  is  the  largest 
eigenvalue  of  the  matrix  M^M.  If  x,  y  €  E°  then  <x,y^>  =  +  ...  + 

xn^n  denotes  scalar  product  of  x  and  y. 

If  f  is  a  function  mapping  a  set  X  into  a  set  Y  then,  for  every 
x  t  X,  f(x)  denotes  that  point  in  Y  into  which  x  is  mapped  by  f;  if  A 
is  a  subset  of  X  then  f[AJ  denotes  the  collection  of  points  in  Y 
iefined  by  f[Aj  =  (y:  y  Y,  y  -■>  f(x)  for  some  x  A). 

We  denote  ty  C(  (-«,0j,  En  ),  or  sometimes  by  C,  the  space  of 
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*  continuous  functions  mapping  the  intervul  (-<»,o]  into  En.  The  topology 

on  C  is  taken  to  be  the  compact  open  topology.1  It  is  fairly  easy  to 

show  (see  Appendix  B)  that  the  topological  space  C  with  the  compact  open 

topology  is  metrixable,  with  metric  p  defined  as  follows:  For  a  fixed 

real  number  b,  0  <  b  <  J,  and  for  a  sequence  of  points  ft^),  0  =  tQ  < 

t.  <  ...  <  t  <  t  <  ...,  with  lim  t,  =  oo,  define,  for  every  q),  ^  in 
L  k  k’»  l  .  k 

_  oo  k-» »  _  _ 

C,  p(<p,  \lr)  =  2  n^,  where  =  minfbK,  supf||cp(t)  -  >|r(  t )  || :  -  t^+1  <  t  < 

k— 0  2 

-t  )  )  C  is  a  complete  linear  metric  space.  It  is  clear  that 
K 

convergence  in  the  compact  open  topology  is  equivalent  to  uniform 


The  compact  open  topology  for  C  is  constructed  as  follows:  For  each 
subset  K  of  ( ~oo, 0 ]  and  each  subset  U  of  En,  let  A(K,U)  denote_the  set 
of  all  members  of  C  which  carry  K  into  U;  that  is,  A(K,U)  =  (qp: 
qp  C  C,  cp[K]  c  u}.  Let  the  family  G  of  all  sets  of  the  form  A(K,U),  for 
K  a  compact  subset  of  (-a>,0]  and  U  open  in  En,  be  a  subbase  for  a  topology 
for  C.  The  topology  for  C  which  is  uniquely  determined  by  this  subbase 
(that  is,  the  smallest  topology  for  C  which  contains  0)  is  called  the 
compact  open  topology.  It  is  denoted  by  3  .  The  family  of  finite 
intersections  of  members  of  G  is  then  a  base  for  the  compact  open 

m 

topology;  each  member  of  this  base  is  of  the  form  H  A(K. ,U.),  where 

i=l  11 

each  K.  is  a  compact  subset  of  (-oo,0]  and  each  U.  is  an  open  subset  of 
E  .  For  a  complete  discussion  of  this  topology  the  reader  is  referred 
to  reference  [21]  and  Chapter  7  of  reference  [22]. 

2 

One  might  hope  that  a  norm  could  be  defined  on  C  so  that  C  would  be  a 
complete  linear  normed  space  (a  Banach  space). _  A  natural  attempt  to 
define  a  norm  might  be  made  as  follows:  Let  ||cp||  =  p(qp,0)  cp  £  C. 

This,  however^  does  not  define  a  norm,  since  the  homogeneity  condition, 
||oscp|]  =  |a|  •  ||cp||  is  not  always  satisfied.  Take,  for  example,  cp(t)  =  0 
for  -oo  <  t  <  -t*T  and  cp(t)  =  (2/t^)t  +  2  for  -t.  <  t  <  0,  and  let  a  =  l/2. 
Then  ||cp||  =  1,  but  ||a  cp|J  =  1  /  l/2.  That  a  norm  cannot  be  defined  for 
C  is  shown  by  Arens  [21]  to  follow  from  the  fact  that  the  domain  of  the 
functions  in  C,  [t:  -oo  <  t  <  0 },  is  not  compact.  The  concept  of  a  linear 
metric  space  is  a  specialization  of  the  concept  of  a  linear  topological 
space.  The  concept  of  a  linear  normed  space  is  a  further  specialization. 
Fortunately,  for  our  purposes  it  will  be  of  no  real  consequence  that 
C(  (-°°,0],  En)  is  not  a  Banach  space.  See  also,  [l8  pp.  ^9-50  and 
PP-  396-397] 
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convergence  on  all  compact  subsets  of  (-®,0]j  in  fact,  a  sequence 

<Pn  -♦  9  in  C  if  and  only  if  for  every  nonnegative  integer  N, 

max{  ||^n(t)  -  <p( t )  || :  -N  <  t  <  0  }  0  as  n  -♦  <».  For  a  given  positive 

constant  H,  we  use  the  notation  C„(  (-®,0],  En  )  to  denote  the  set 

{9:  9  £  c>  sup{  ||<p(  t )  |( :  -«  <  t  <  0  )  <  H  J.  Again,  we  shall  abbreviate 

the  notation  to  Cjj  when  the  meaning  is  clear. 

Let  A  >  -  00,  and  let  x  be  a  continuous  function  mapping  the  interval 
(««,A)  into  En.  Then,  for  every  t,  -«  <  t  <  A,  we  denote  by  x^  the 
translation  to  the  interval  (-»,0],  of  the  restriction  of  x  to  the 
interval  (-«,tj;  that  is,  x^  is  an  element  of  C,  defined  by  x^(a)  = 
x(t  +  a)  for  -00  <  o  <  0.  In  other  words,  the  graph  of  x^  is  the -graph 
of  x  on  (-®,t]  shifted  to  the  interval  (-®,0]. 

If  r  is  a  real  number,  if  f  is  a  function  mapping  CH  into  En,  and 

if  x(t)  denotes  the  derivative  of  x  at  t  >  r,  we  consider  the  following 

autonomous  functional -differential  equation: 

x(t)  ?(xt),  t  >  r  ,  (2-1) 

Equation  (2-1)  is  called  a  functional-differential  equation  because 
each  element  of  the  vector  x(t)  is  determined  by  the  value  of  a  functional 
on  C.,.  We  say  that  x(r,m)  is  a  solution  of  Equation  (2-1)  with  initial 
.■condition  m  t  C„  at  t  -  r  if  there  exists  some  A  >  r  such  that  x(r,<p) 

Is  a  mapping  from  (-°°,A)  into  En  with  xt(r,9)  in  C}1  for  r  <  t  <  A, 

7  ,'r,5)  -  9,  and  if  x(r,cp)  satisfies  Equation  (2-1)  for  r  <  t  <  A. 

The  concept  of  a  functional-differential  equation  is  more  general 
■  .an  that  of  an  ordinary  differential  equation.  Consider  any  ordinary 
inferential  equation,  for  example  x(t)  =  2x(t)  +  x^(t).  If  xi,t)  is  a 
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solution  of  this  equation,  the  derivative  of  x  at  some  point  t  may  be 
computed  if  one  knows  only  the  value  of  x  at  the  point  t.  For  a 
functional -differential  equation,  the  value  of  the  derivative  of  a 
solution  at  some  point  t  depends  upon  the  values  that  the  function 
x  assumes  over  an  interval  for  which  the  point  t  is  the  right-hand 
end.  Clearly,  the  concept  of  a  differential-difference  equation  is 
also  a  special  case  of  a  functional -differential  equation.  From  this 
point  of  view  it  is  obvious  that  an  initial  condition  for  a  functional- 
differential  equation  and  also  the  state  at  some  time  t  of  a  system 
which  is  governed  by  a  functional-differential  equation  should  be 
specified  by  a  point  in  some  space  of  functions. 

In  a  manner  similar  to  that  used  for  ordinary  differential  equations 
and  differential -difference  equations  (see  references  [23]  and  [2h]), 
one  may  prove  the  following  existence  and  uniqueness  theorem:  If  f 
is  continuous  in  C^,  then  for  any  qp  in  there  is  a  solution  of 
Equation  (2-1)  with  initial  condition  qp  at  t  =  r.  If  f  is  locally 
Lipschitzian  on  CH;  that  is,  if  for  any  Hx  <  H,  there  exists  a  constant 
L(HX)  such  that  ||f(<p)  -  f(f)||  <  L(H1)p(qp,^)  for  all  <p,*  in  CH  with 
p(i,0)  <  p(*,0)  <  H1,  then  there  is  only  one  solution  with  initial 

condition  <p  at  t  r  and  the  solution  x(r,<p)  depends  continuously  upon 
<p.  Also,  f(cp)  locally  Lipschitzian  in  qp  implies  that  the  solution  can 
be  extended  in  C  until  the  boundary  of  is  reached. 

A  set  M  in  C  is  called  an  invariant  set  if,  for  any  <p  £  M  there 

exists  a  function  x  defined  on  (-»,oo)  with  (L  M  for  every  t  in  (-«o,«o) 

and  xq  -•  <p,  such  that,  for  every  o  in  (-«,»),  if  x  (o,x^)  is  the 

solution  of  Equation  (2-1)  with  initial  condition  x^  at  o,  then 
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-* 

‘  xt  "  Xt  f0r  ^  >  O' 

If  V  is  a  continuous  functional  on  CH,  and  if  x(0,q?)  is  the 
unique  solution  of  Equation  (2-1)  with  initial  condition  9  at  t  »  0, 
we  define  V^g^Cqj)  and  V*2-l)^by: 

V(2_i)(9)  »  lira  +  £(  V(xh(0,9)  "  V(9)  ), 
h-*  o  • 

^(2-i)(^)  3  —  ■  v(<?)  )• 

h-»  o+ 

2.  Stability  Theory  for  Functional -Differential  Equations 

If  f(0)  a  0,  then  the  solution  x  **  0  of  Equation  (2-1)  is  said  to 

be  stable  if  for  every  e  >  0  there  exists  a  6  >  0  such  that  <p  G  CH  and 

p(q>,«)  <  6  implies  that  x  (0,<p)  exists  for  all  t  >  0,  is  in  C„,  and 
p(xt(0,qj),  0)  <  C  for  all  t  ^  0.  If,  in  addition,  there  exists  a  5  >  0 
such  that  p(qj,0)  <  6  implies  that  xt(o,9)  is  in  CR  for  all  t  >  0  and 

x,.(0,q))  -*  0  as  t  -♦«,  then  the  solution  x  =  6  is  said  to  be  asymp- 

totically  stable.  If  the  solution  x  =  0  is  asymptotically  stable  for 
all  H  >  0  and  all  5  >  0,  then  the  solution  x  =  6  is  said  to 

be  completely  stable.  The  solution  x  =  0  is  said  to  be  unstable  if  it 
is  not  stable. 

It  should  be  noted  that  if  x  is  a  continuous  function  from  (-»,<») 
to  En  then  lim  ||x(t)||  =  0  if  and  only  if  lim  p(xt,0)  =  0  since  con- 

t<-*  <0  t-4X» 

vergence  in  the  compact  open  topology  is  equivalent  to  uniform  con¬ 
vergence  on  all  compact  subsets  of  (-®,0],  in  particular  the  set  (0). 
Thus,  defining  stability,  etc.,  in  terms  of  the  compact  open  topology 
yields  the  desired  properties. 
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We  now  state  three  theorems  due  to  J.  K.  Hale  [12]  concerning  the 
stability  of  equilibrium  solutions  of  functional-differential  equations. 

As  Hale  points  out,  these  theorems  generalize  the  results  of  LaSalle 
[25,263  for  ordinary  differential  equations.  The  proofs  are  also  extensions 
of  the  ones  given  by  LaSalle  and  are  to  be  found  in  Hale’s  paper.  The 
proofs  given  by  Hale  in  [121  are  actually  stated  for  functional-differential 
equations  on  the  space  C(  [-r,0],  En),  where  r  is  a  positive  real  number; 
however,  as  is  pointed  out  in  the  last  section  of  his  paper,  the  theory 
applies  as  well  to  functional-differential  equations  on  C(  (-00,0],  En). 

In  some  places  the  wording  of  these  theorems  has  been  changed  slightly; 
and  in  Theorem  3  (Theorem  1*  of  Hale's  paper)  a  trivial  change  has  been 
made  in  condition  l) .  In  all  of  the  following  theorems  we  assume  that 
the  function  f  in  Equation  (2-1)  is  continuous  and  Lipschitzian  on  Cjj, 
for  H  >  0. 


Theorem  1.  Let  C*  =  U  C  .  Let  V  be  a  continuous  functional  on  C*.  If 

0<r<o°  r 

designates  the  region  of  C*  where  V(cp)  <  £,  suppose  there  exists  a  non¬ 
negative  constant  K  5  ||cp (0 )  |l  <  K,  V(cp)  >  0,  and  V^_1j(cp)  <  0  for  all 
cp  €  U  .  If  R  is  the  set  of  all  points  in  where  V^-l)^)  =  ®  an<*  ^ 
is  the  largest  invariant  set  in  R,  then  every  solution  of  Equation  (2-1) 
with  initial  condition  in  U  approaches  M  as  t-»  00. 

Jv 


Theorem  2.  Suppose  f(0)  =  0,  and  let  the  continuous  functional  V  be 

defined  on  C*  =  U  C  such  that  V(o)  =  0.  Let  U  denote  that  region 
_  o<y<M r  1  _ 

of  C*  where  V(cp)  <  £.  Assume  that  there  exists  K  such  that  ||cp(0 ) ||  <  K 

for  all  9  €  U  .  Let  u(s)  be  a  function,  continuous  and  increasing  on 

a 

[0,K ),  where  u(0)  =  0.  If  CH  <=  u  ,  if  u(|ji(0)||)  <  V(i),  and  if  V^2-1j(i)  <  0 
for  ’all  cp  €  U  ,  then  the  solution  x  =  0  of  Equation  (2-1)  is  stable. 

a 

15 


Theorem  3-  Suppose  ? (0)  ■  6  and  V  is  a  bounded  continuous  functional1 
on  d  and  there  exists  a  y  and  an  open  set  U  in  C  such  that  the  following 
conditions  arc  satisfied: 

1)  V($)  >  0  on  U  n  C  ,  V{£)  u  0  on  that  part  oi  the  boundary  of 
U  in  CYi 

2)  6  belongs  to  the  closure  of  U  fi  C^i 

3)  v'$)  <  u(  15(0)1!  )  on  U  fl  C  ;  where  u(s)  is  continuous,  non- 
negative  and  nondecreasing  on  (0,H),  and  u(0)  *  0> 

k )  V(2  5  0  on  the  closure  of  U  A  and  the  set  R  of 

•  *  .  _  , 

<p  in  the  closure  of  U  n  such  that  ^ j ( qp )  =  0  contains  no  invar¬ 

iant  set  of  Equation  (2-1)  except  <p  =  6. 

Under  these  conditions,  the  solution  x  =  0  of  Equation  (2-1)  is 
unstable  and  the  trajectory  of  each  solution  of  Equation  (2-1)  with 
initial  condition  in  U  n  C  must  leave  C  in  some  finite  time. 

r  r 

^By  bounded  we  mean  here  that  tnere  exists  i  >  0  3  V(<p)  <  ■£ 

for  all  ^  c  Cjj- 
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TRANSMISSION  LINK  FUNCTIONAL  REPRESENTATION 

In  this  chapter  we  shall  develop  several  pairs  of  functional 
equations  which  describe  the  electrical  behavior  of  certain  simple 
two-port  networks  which  contain  a  lossless  transmission  line.  The 
two-port  networks  which  we  shall  be  concerned  with  are  constructed 
by  connecting  a  resistor,  either  in  series  or  in  parallel,  at  each 
end  of  a  lossless  transmission  line.  We  consider  such  two-port 
networks,  rather  than  simply  treating  the  transmission  line  itself 
as  a  two-port  because,  if  it  is  required  that  at  least  one  of  the 
resistors  have  a  finite  value,  greater  than  zero  (which  shall  be 
the  case  in  tne  app  ication  of  these  results),  then  the  linear 
functionals  which  occur  in  our  resulting  equations  will  be  contin¬ 
uous.  That,  the  functionals  have  this  property,  in  fact,  that  the 
functionals  satisfy  a  Lipschitz  condition,  will  be  proved  in  the 
final  section  of  this  charter. 

.1  ■  The  Transmission  Line 

We  define  a  lossless  transmission  line  to  be  a  distributed 
electrical,  two -port  network,  as  shown  in  Figure  3-1,  which  is  charac¬ 
terized  by  three  parameters:  l,  the  length  of  the  transmission  line; 

C.  the  distribute  i  capacitance  per  unit  length;  and  L,  the  distributed 
inductance  per  unit  length  We  always  assume  / ,C ,1  >0.  It  is  con¬ 
venient  to  define  twe  additional  parameters  which  may  be  used,  along 
wit;.  I ,  to  give  an  alternate  method  of  characterizing  the  line.  We 
lefine 
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Figure  3.1.  A  typical  lossless  transmission  line. 


ao  =  V  L/c  , 

a  «  V  in'  . 

ZQ  is  called  the  characteristic  impedance  of  the  transmission  line  and 

a  is  the  reciprocal  of  the  line's  phase  velocity.  We  also  define  Yq, 

the  line's  characteristic  admittance,  by  Y  =  1/Z  . 

o  o 

If  x  is  a  point  in  the  interval  [0,i]  and  t  is  a  point  in  the 
interval  (~*>,Tj,  for  some  T  >  then  at  any  time  t  and  at  any  point 
x  the  transmission  line's  voltage  and  current  are  denoted  by  v(x,t) 
and  i(x,t)  respectively  Thus,  v  and  i  are  real  valued  functions 
defined  on  the  set 

Hj,  --  (  (x,t)  :  0  <  x  <  X,  -ae  <  t  <  T  j 
I’ote  that  v  (t)  »  v(o,t),  i  (t)  -  i(o,t),  v  (t)  =  v(f,t),  and  i  ( t )  = 
-i(X,t).  We  assume  that  for  each  lossless  transmission  line  in  our 
distributed  networks  there  exist,  functions  f^  and  such  that 

v(x,t)  f^(ax  -  t)  +  f0(ax  +  t), 
i(x,')  ~  [  fj (ax  -  t)  -  f  (ax  +  t)  J, 
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{3-1} 


for  *11  (x,t)  €  My.  From  these  equations  we  obtain; 
v(0,t)  -  f^-t)  ♦  f£(t), 


Ko, t)  =  fx(-t)  - 1~  f2(t), 

(3-2) 

O  0 

v(i,ai-t)  =  fL(t)  +  f2(2al-t), 

(3-3) 

1(1, al-t)  =  fx(t)  -  f2(2al-t). 

<3-M 

From  Equations  (3-3)  and  (3-M  we  find 

fx(t)  -  |[  v(l,  ai-t)  +  ZQi(i,  ai-t)  ],  ( 3-5) 

and  from  Equations  (3-1)  and  (3-2)  we  have 

f2(t)  =  §[  v(0,t)  -  ZQi(0,t)  ).  (3-6) 

Substituting  Equations  (3-5)  and  (3-6)  into  Equation  (3-M  and 
replacing  t  by  2ai-t,  then  substituting  Equations  (3-5)  and  (3-6)  into 
Equation  (3-1)  and  replacing  t  by  t-ai  gives 

v(0,t)  -  ZQi(0,t)  »  v(i,t-ai)  -  Zui(i,t-a i) 

and 

v(i,t)  +  ZQi(l,t)  -  v(0,t-a i)  +  ZQi(0,t-a/). 

Thus,  letting  t  =  a t  and  recalling  that  i(i,t)  =  -i#(t), 

£ 

vo(t)  -  ZQi0(t)  =  v^(t-r)  +  Z^U-r), 

v/t)  '  V/(t)  =  V0(t_T)  +  Vo(t^)'  (3-7) 

Equations  (3-7)  are  used  as  the  starting  point  in  the  next 
section  where  the  functional  equations  which  govern  the  behavior  of 
the  various  two-ports  are  derived. 
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2  •  The  Transmission  Line  Two -Porta 

In  this  section  Equation*;  (3-7)  are  used  to  derive  the  functional 
equations  which  describe  the  electrical  behavior  at  the  ports  of  the 
three  two -port  networks  shown  in  Figure  3.2.  That  is#  for  example. 


0^  0 — ^ 

v  <  V 

°(  0 

Z0,a,| 

\'l 

V 

A .  . . 

J 

>  ) 

(a) 

Figure  3.2.  Transmission  line  two-ports 


for  any  t  6.  (-»,T]>  T  >  we  will  show,  for  the  network  of  Figure  3. 2a 
how  i  (t)  and  i  (t)  may  be  expressed  as  functionals  whose  arguments 
are  the  functions  (  of  0,  -»  <  0  <0),  v  (t+o)  and  v  (t+ci).  The 

“  O  It 

resulting  functional  equations  are: 
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For  Figure  3.2b: 
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For  Figure  }.2c: 
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;  *  ( 2k+i  )t ) » 

t)k  v/(t  -  2(k+l)t) 
t  -  (2k+l)r)  . 

{)k  vo(t  -  2(k+l)t. 

■  (2k+l)T  , 

i^(t  -  2( k+1 )x | 

-  (2k+l)t/. 

K  iG(t  -  2(k+l)t 
-  (2k+i),  , 
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v,<t> "  1,U)  * 2 


Vo 


<«,♦»<,>  k» 


X  {tWk  ‘  2<k+i*) 


*  s  K\  tJEo.  7  ,  )  £  (V/  ‘o'*  •  < ’•«> 

kaO 


In  Equations  (3-6)  Tq  and  are  the  usual  reflection  coefficients 
at  each  end  of  the  transmission  lines  when  the  independent  port 
variables  are  set  equal  to  zero;  e.g.,  for  Figure  3.2a,  VQ  » 

(H0  -  Z0)/(R0  +  ZQ)  and  »  (R^  -  Z0)/(R|  +  ZQ).  We  assume  for  each 
network  that  both  resistors  are  nonnegative  and  at  least  one  has  a 
positive  real  value  (hence  |rQr  J  <  l). 

Equations  (3-7 )>  of  course,  apply  to  the  primed  variables  in 
Figure  3 -2a*  The  primed  variables,  however,  are  related  to  the 
unprimed  variables  by  the  following  equations: 


^  -  V 


V  »  v  -  R  i  , 
O  O  O  o' 


l,  .  1,, 


yj  '  vi  •  Vi- 


Thus,  Equations  (3-7)  may  be  written  in  terms  of  the  unprimed  variables 


as 


Vl)  •  Vo(t)  -  V»(l>  '  V-y>  ■  V,(t-,)  * 

•  V,U)  -  V(l,!  V-'i  •  V.'1-’1  *  zoVt'r>' 

u  ni  hence , 

'  "Tvl  V ’ "  »7s  '*U'T)  *  Ro  *  1 
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1 


( 5-9) 


R  m  Z 

i/t)  -  ■  v  (tj  -  -  vo(t-t)  ♦  i0(t-r). 

Replacing  t  by  t-r  in  Equations  (3-9)i  have 


R ,  -  SL 


■  irh;  -  irii;  v,<t-£*>  *  sf-rs* 

‘i1*-13  ■  St4t  Vt',)  -  B-iT  vo(t"8«>  *  r-r*2  V*-8*)-  <>-“3 

4  0  4o  to 

Substituting  Equations  (3-10)  into  Equations  (5-9)  givt/a 

1o(t)  “  rr\'  vo(t)  *  r«7^o^i  v*o)  v0(t"2t) 

H.  -  z_  T  (R,  -  Eo)(Ro  -  Zo} 


R  +  Z 
o  o  *■ 


.  -jL 


*/  +  Z0  i 


v4't“T)  +  (fTTz^  ♦  z"7  io(t*3t)' 


(R  -  2.) 


15  i^h-  v.(t)  •  TsjT^n^-^i  v/t*2T) 


r„  +  z„ 

i  0 


1  - 


F  -  Z 
_o  _o 

H  Z„ 
o  c 


(Rc  -  z  }(R  -  Z  ) 

V*-'  *  tSTTz^Tz^y  V1**3' 


and  hence, 


‘o'1' '  «rH  V-'3  •  rrr  r<  V1-2"3 

C  O  OO 

■  2(R0  +  Zo)(Pz  +  ZC}  *  r°ri  io(t_2T)' 

R,  +  Y  Vf^  ‘  R  +  z  r0V£('t'2x^ 


-  ?  rB—/iWprro)  vo(t-*>  +  ror/  V*-2^- 


(3-11) 


We  nov  note  that  the  second  of  Equations  (3-11)  may  be  obtained 
from  the  first  by  simply  excnanging  the  o  and  i  subscripts  on  the  symbols 
i,  v,  R  and  I’.  Due  to  this  symmetry  ve  need  nov  consider  only  the  first 
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of  If  r  /  C  «md  V ^  $■  0  tActl,  tQT  k  “ 

\k  ,  / .  n(,  *  d 

’o  *  -o 


»VV‘  ‘o''-21"’  -  ox  eyy*  t,"-2"’  -  rrr0 


■2  . ';rjt»rr®  (tw  v»(t-<2Mi>,>  *  (iw 

therefore ,  for  p  =*  1  >  2 ,  .  •  •  < 

p  ,  p 

;  <|v'<>k  ‘0"-2k'>  -  ox 

k*o 

P 


k+1 


Y  <'Wk  ‘o"-2k,)  -rr-r  t  <ror/  ,o('-2k,) 

K  -o  °  k~o 


fr  Z  W1  vo(t-2(k+1)t) 


o  o 


k*o 


-  X  «  zife;  rr;  Z  <v/  V‘-'2kkl>k> 


o  o'*  /  o 


k-o 


£  (r0rjf)k+i  io(t-2(k+i)t). 


k^o 

But  i 

P 


P*1 


jv/*1  ‘0(p-pU«),)  ,  JVr.r/  i^t-pk,); 


o 


k=l 


!.«r«I'ure , 

>  T Z|r«r/  vo(t- 

o  o  ^ 


k-1 


k-o 


TY"  £«  1  v0U-2(k+ib)  -  2(f  +  Zo)(Rf-^: 

3  o  * — ■  0  0  10 

k-o 


r\ 


7T  ♦  uy,)p+1  i0(t-^(p+i)T). 

2k 


0,  1,  2,  ...  , 

l0(t-2(k+l)t)i 


2kT ) 

P 

iJ„yVkv/(t.(,k.!), 

k--o 


Also 


r  s 

DW  v0(t“2kt)  »  v0(t)  ♦  £  (rol|)k+i  vo(t.2<k+l)r) 


-  ^cri)p+i  v  (t.2(p,i}t), 


therefore , 


>(t)  *R7T-2“[Vo(t5  +  PC1  *  ro>  ^  V0(t-2(KH}t)] 

k«o  ^ 

zo  P 

'2  <"0  +  Z0)(Ri  +  V  ^  (WK 


(r3r/)p+1  V0(t-2(P+1)T)  ♦  (rQri)p+1  iQ( t-2(P+i. )T ) . 


Using 


r£+1  -  ro  *  r£  <r0  *  D  *  -2  rk  , 


0  o 


we  obtain 
i0(t)  aRl 


XX  h"1 " 2  "T^X  I/W* 

k«o 

zo  P 

2(R0  +  V(Ri  +  V  Pv/ 

k-o 


"  R0  +  Z0  ^r°rjl^P  1  V^<P+^>  +  (WP+1  i0(t‘*-J(P+1)T)  ■ 

Now,  since  0  <  R  <  ®  and/or  0  <  H.  <  -,  and  hence  jr  r.l  <  ]  . 

*  o  l 

we  have,  if  vq,  iQ,  v^,  are  bounded  on  ( -oo ,T 3  • 

n"  i~h  Cor  J"*1  v  (t-2(p.i)T)  o, 

p-«  a,  0  O 
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lifl  loU-2(i*l)*)  -  °, 


r 

U»  v0(t-2(**J-)‘r)  “  £<V/  v^t-^k+ljt),  etc ., 


converge-  Hence ,  using  the  symmetry  referred  to  following  Equations 

(  Vll),  we  have  derived  the  first  two  of  Equations  (3-8).  Immediately 

following  Equations  (M)  it  was  assumed  that  I’Q  4  0  and  1'^  ^  0. 

We  may,  however,  compare  Equations  0-11}  with  the  first  two  of 

Equations  (3-8)  and  see  that  the  latter  are  valid  for  all 

I'c .  t‘x  (  |rQ|  <1,  |l'x |  <1,  |rQrx|  <  1  ),  provided  we  define 

(r  r,)°  »  J  if  r  =  C  or  r.  =  o.  We  shall  make  this  definition, 
o  Jt  o  i 

The  remaining  four  of  Equations  0-8}  are  easily  derived  in  a 
•  similar  manner. 


y.  The  Lipschitz  Condition 

In  this  section  we  prove  that  the  functional  equations  describing 
the  electrical  behavior  at  the  ports  of  the  networks  considered  in 
'he  last  section  satisfy  a  Lip:; -'hit 2  condition  provided  that  the 
compact  open  topology  on  C  is  metrlzed  with  an  appropriate  metric. 
a’c  consider  the  mapping  f,  from  C’  (  (-®,0j,  t'n  )  into  E  ,  defined 
as  follows:  For  i  -  1,  ...  ,  n,  let  >0,  rJ1^  >0,  1  ^  >0, 

Z  S  '  >  0,  and  let  t ' 1  ^  ^  f  1  \  +  3* 1 '  for  k  »  1 ,  s  , 

A  /.-l  K  K-a 

Let  A  be  a  real  i.-t  diagonal  matrix, 
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Where  jaj 

A 

►— 

c 

a.  f  *  *  *  $  n  *  Lc  1  M|  1 

M.)(  M,  denote  real  n 

-  J 

matrices . 

Then,  w  9  < 

V 

»  m 

?(q)  a 

^9(0)  *  h2 

l  %  * «,  I 

.k  - 

A  *  , 

(3*12) 

k=o  k=o 

•V 

where  o 

denotes  the 

vector 

_  (  t(n)l 

»  •  *  *  '  V  \  > 

J1  , 

and  similarly  for  $>a  .  If  we  consider  any  finite  collection  of  two- 
port  networks  of  the  type  considered  in  the  last  section,  the  system 
of  functional  equations  describing  the  electrical  behavior  at  their 
ports  is  of  the  form  specified  in  Equation  (3-12)  above. 

We  first  define  the  metric  p  with  which  the  compact  open  topology 
on  C  will  be  metrized:  Let  a  =  max  {  ^ j  :  i  =  1,  ....  n  ),  and 

choose  b30<a<b<l  Let  tQ  «  0 ,  t^  >  max  { O  :  i  <  1 ,  .  .  •  ,n  J , 

*  <11(11  .  #  #  « 
and  pick  T  >  max{  T'  ,  S'  .  i  *  1,  •  •  •  ,  n).  Let  t^  =  *  T 

for  k  i  2,  3,  ...  .  For  i  =1,  ...  ,  n,  define  the  integers1 

=■  [iV11]  ♦  1,  "  »  £t#/S^|  +  1.  Clearly,  each  interval 

*  * 

(-t^^,  -t^J,  for  k  '<  1.  2,  ...  ,  contains  at  least  one  of  the  points 

and  at  moat  N '  such  points  (for  1  1 ,  ...,  n);  and,  similarly, 

each  such  interval  contains  at  least  one  of  the  points  -s^1 ^  and  at 

most  ll ^  '  such  points  .  If  * "  -  [  ( t^  -  ^  )/  1 '  J  ♦  1 ,  and 

f,i  * ' "  "  (^l  *  ♦  1,  th"n  the  interval  (-t*,  Oj 

contains  "  of  the  points  -t^^  ant*  of  the  points  -sj^  , 


4 We  use  here  the  following  notation:  If  r  is  any  real 
number  then  ,r.  denotes  the  greatest  integer  k  sue;.  ?*  ••»  k  <  r 
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for  i  «  1,  . . . ,  n.  Let  N  ■  max  (  ,  H'' ,  "#»4 

*  * 

then,  each  interval  [-t^,  -t^J,  for  K  3  0,  1,  S, 
least  one  of  the  points  -tj^  ^  and  -sj^  and  iot 


:  i*l,  ...»  n J ; 
}  contains  at 
'  N  such  points, 


for  i  =1,  •••>  n.  Let  the  compact  open  top  v  on  C  be  metrized 


with  thf  metric  p  defined  by  p(^,j)  =  2  m.  ,  wi_  e  bl» 

K=o 

jrin(  bk,  sup {  |^{t)  -  *(t)||  :  -t*+1  <  t  <  -t£)  },  V  <£,  f  €  C. 

We  now  prove  two  lemmas  which  will  be  used  in  our  procedure  for 


specifying  the  Lipschitz  constant  for  f . 


Lemma  1 .  if  -t '  denotes  an  arbitrary  fixed  point  in  one  of  the 

#  * 

intervals  (-t^+1  ,  -t^]  and  if  denotes  the  mapping  from 
CH(  (-**,0],  En  )  to  En  defined  by  g^p)  *  $(-t^),  V  i  c  Cy,  then  ^ 
satisfies  a  Lipschitz  condition  on  with  Lipschitz  constant  = 
max(  1,  2Hb“k  ). 

•  _  ,  ■  j. 

proof:  Let  q>,  f  €  C^-  If  ||<p( -t^)  -  *(-tk)|l  <  b  ,  then 

li  q(-\)  -  t(-tj[)||  <  <  1^  p($,f).  It  l^(-t^)  -  *(-\)|; 

>  b‘l  then  p($,f)  >  m^  =  bk.  But  then,  b_kp(qi,*)  >  1  =>  2Hb‘kp(9,f )  > 

2H;  or  (since  |l  ip(  -t^)  -  ♦(  -tR) ||  <  2H),  ||  q>(-tk)  ~  *(-tk)il  < 

2Hb"k  p(<jj,f)  <  L  p($,5).  Q.E.D. 

Lemma  2.  If,  for  k  a  0,  1,  2,  . K,  -t^  denotes  an  arbitrary  fixed 

*  # 

point  in  the  interval  [-tk+1>  -t^j,  0^  denotes  a  real  number,  and  g^ 

denotes  a  mapping  of  C„(  (-*,0J,  En  )  into  En  as  described  in  Lemma  i ,  then 
K  a 

if  Ly  =  2  i^l'^k'  where  is  the  Lipschitz  constant  for  each 
k— 0 

(which  is  guaranteed  to  exist,  by  Lemma  1),  then 
K 

o(^i  +  )>  *  ?,♦  £  CH. 
k-  o 
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*Mat!  £  .  w-t,-)!  <  KHiii-t,;)  -  s 

K 

i  Ki’1^  =  Sc  P(<P>i')>  V-  c  CH.  q.E.C. 

It.  should  be  noted  that  the  constants  1^  and  of  Lamas  1  and  2 

do  not  depend  on  the  particular  point  -t^  specified  in  the  interval 

t-tk+l'  'V' 

We  now  specify  the  Lipschitz  constant  L  for  the  mapping  ?;  that 
is,  we  determine  the  constant  L  such  that  #*($)  -  f(f)||  <  L  y 

<p,fe  CH:  Pick  the  positive  integer  J  such  that  J  >  J=»(a/b)^  < 

2HNn  and  <  "ito  *  Let  LJ-1  denote  the  constant  specified  in 
Lena*  2  for  the  sequence  of  mappings  (g^  and  constants  =  a^n, 

k  "  0'  J"1’  and  let  *0  deoote  the  Upschitz  constant  specified 

In  Learn  1  for  gQ<<p)  =  <j5(o).  Then,  let  L  *  {(Mj^ j]  - ^  + 

(  11*211  *  IN,»Kl  +  Lj^). 

Theorem.  If  L  is  specified  as  above  then  ¥  $,f  £  j|?($)  -  ?(*)| j  < 

L  p(<P,f).  vhere  f  is  the  mapping  from  CH(  <— ,0],  En  )  to  En  defined 
by  Equation  (3-12). 


Uroof:  ¥  <p,*  e  cH: 

n  qq 

l'f('P)-  ?(*)ll  <  11^(0)  -  ^(0)11  +  £ 

k=o 


.k,  - 

A  (<p. 


k-o 


*t  )! 

k 


+  i:m. 


7 

t _ _ 

k=c 


Ak(i 


-  f 


»!  <  l.!MLl!-iop(<p,t) 

K 
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w-«  '  *Kn  f  -  *s  )#* 

‘  x  k  k  k 


We  now  show  that  each  of  the  terms  }|^  AK($.  -  f.  )|j,j|  W<®  -  *  )j! 

\  \  L*  Bk  Bk 

kao  *  k=o  k  k 

Is  le  <»  than  or  equal  to  (l  +  Lj _^)  p(<p,}),  and  then  the  theorem  is 


proved- 


t  A‘<\  -  »t  >«  s  £  -  M  -  f>k 


-  *t  II  -  K.  -  v  II  < 

“  *  kaO  k  k 


t/{  )i )  <  Z  .“(  -  w-4‘> 


k-O  lal 


t  ( i  ■kw-4°)  -  «- 


tij)ii 


1=1  k=0 


For  each  i  -  1,  ...  ,  n,  3  a  subsequence  [-t^^J,  j  =  0,  1,  2,  ... 

J 

consisting  of  exactly  one  point  from  each  interval  ,  -t*| 

at  which  ||<p(-t^i ^ )  -  ^(-t^^)jj  is  u  maximum  for  all  of  the  (at  most  N) 
points  -t^  in  the  interval.  Obviously, 

CD  015 

^ak|'9(-t[i))  -  *(-t[i))||5  -  *(-tj‘i))||  for  i=l, . . .  ,n.  TV. 


llZ  A*{\  -  \)n  <  Z  ( Z  - 

k-o  *  k  i=l  ‘  j=o  J 

j^O  *i=i  J  J 


Tor  each  J  *  0,  1*  2,  *• 

* 

in  each  interval  [-tfc+1. 


,  3  one  point*  of  the  n  points  {-tj^si-l, . 
■t*J,  at  which  || )  -  f(-tii^  )||  is  a 


>«} 


n 

maximum.  Call  it  -tj.  Then,  £  |ft(  -  $( -t£l  * )  |  < 

i-1  ^  J 


aM-tj)  -  *(-tj)||  and  hence, ||£  Ak($t  -  ^  )||  <  y*aJlln||(p( -tj)  -  f(-tj)|| 

k*>  k  k  J-o 

J-l  • 

-=  £  ^Nn'I^-tj)  -  J(-tJ)|j  +  £  aJNn|ft<-tj)  -  f(-tj)||  < 

J=o  J=J 

«D 

+  a^Nn||^(-tj)  -  J(-tj)||,  by  Leona  2.  But,  since  J  >  j  •=>  > 

2HNna^,  a^NnJI^C-t')  -  J(-t')||  <  a^Sn-2H  <  b*^  for  all  J  >  J.  Therefore 
(since  also  a"Hn  <  1),  J  >  J  *=> =  min{b^,  sup(jjq>(t)  -  $(t)||: 

-t*^  <  t  <  -tj)  )  >  min  (bJ,  sup(aJlto|l9(t)  -  ♦(t)||:-t*+1  <  t  <  -t*}  } 

>  min  (a^Nn^-t')  -  f(-t')||,  sup(  aJNn||<p(t)  -  f(t)||:  -t*+1  <  t  <  -t*)  ) 

i  *o  oo 

«  aJNh||$-t')- i(-t')||.  Thus,  p($,t)  -  Z  m  >  Z  m, 

OB 

>  £  aJNn|fc(-tp  -  ♦(-tJ')||. 


Thus , 

00 

|l  ^Ak(9t  -  )||<  Lj^  p(qi,f)  +  p($,^)  =  ^Lj.i  +  p(0>*)' 

k=o  k  k 

The  same  technique  exactly  proves  that 
00 

li  V  AK(<ps  -  5s  )||  <  (Lj  a  +  l)p(ip,f)- 
k  k 
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q.E.D. 


Chapter  IV 


STATE  SjJJATXOaS  FOR  A  CUSS 
OF  NONLINEAR  DISTRIBUTED  NETWORKS 

In  this  chapter  we  define  the  class  of  nonlinear  distributed  network 
for  which  our  stability  criteria  apply.  We  also  show  that  the  behavior 
of  any  network  in  this  class  is  determined  by  a  system  of  functional- 
differential  equations  having  the  form 

x(t)  =  ?(xt),  t  >  0, 

where  the  state  of  the  system  at  time  t  >  0  is  represented  by  x  ,  a  point 

w 

in  the  space  Cjj(  (<*>,0),  E*1  ),  and  the  function  f  is  continuous  and 
locally  Lipscbltzian.  The  form  of  the  function  f  is  also  given,  so 
that  if  one  selects  any  network  from  the  given  class  it  will  be  evident 
how  to  construct  the  particular  functional-differential  equation  which 
determines  its  behavior.  We  first  consider  the  writing  of  state 
equations  for  lumped  networks. 

] .  Lumped  Networks  and  State  Variables 

Much  has  been  written  cn  the  subject  of  writing  state  equations  for 
lumped  networks.  See,  for  example,  references  [27]  through  [36]. 

Although  both  linear  and  nonlinear  networks  have  been  considered,  we 
restrict  our  attention,  for  the  moment,  to  lumped  linear  networks. 

Let  the  integer  n  denote  the  number  of  Independent  voltage  and  current 
sources  in  a  given  linear  network.  We  may  then  consider  the  network 
to  be  a  lumped  linear  n-port,  containing  no  independent  sources,  with 
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independent  voltage  and  current  sourcea  connected  at  each  port.  For 
nofit  auch  networks  it  is  possible  to  designate  a  cert#4n  '-'•'llcction  of 
the  a- port 'a  branch  voltagea  and  currents  es  the  "atate  variables"  of 
the  network.  These  state  variables  have  the  property  that  the  voltage 
across  or  the  current  through  any  branch  of  the  network  has  a  unique 
representation  as  a  linear  combination  of  the  state  variablea  and  the 
independent  source  voltages  and  currents.  Thus,  the  behavior  of  the 
network  is  completely  determined  if  the  behavior  of  the  state  variables 
and  the  Independent  sources  is  known. 

Usually  it  is  possible,  and  convenient,  to  select  as  state  variablea 
the  voltages  across  capacitors  and  the  currents  through  inductors  in  the 
n-port.  We  shall  not  dwell  on  the  matter  of  when  it  is  possible  to 
select  such  a  set  of  state  variables  to  characterize  a  given  linear 
network  since  this  matter  has  received  auch  attention  in  the  literature 
[27,  28,  29,  30,  31,  32]. 

When  it  is  possible  to  characterize  a  linear  network  having  no 
mutually  coupled  inductors  by  a  set  of  state  variables  as  described  above, 
we  may  write  the  system  of  linear  differential  equations 

P  x(t)  -  A  x(t)  +  0  u(t),  t  >  0,  (4-1) 

where:  x(t)  denotes  the  state  vector,  a  vector  whose  components  are  the 
state  variables  for  the  network  (the  voltages  across  capacitors  and  the 
currents  through  inductors);  u(t)  denotes  the  vector  whose  components 
are  the  values  of  the  independent  sources;  P  denotes  a  diagonal  matrix 
in  which  each  diagonal  element  is  equal  to  the  value  of  the  reactive 
element  (capacitor  or  inductor)  corresponding  to  the  i-th  state  variable. 
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We  note  that  P  Is  a  positive  definite  sytmetrie  matrix.  Thus,  the 
left-hand  side  of  Equations  (k-l)  is  equivalent  to  a  vector  whose 
components  are  equal  to  the  values  of  the  currents  through  capacitive 
branches  and  the  voltages  across  inductive  branches  in  the  network. 

The  right-hand  aide  of  Equations  (4-1)  is  the  expression  of  these  sane 
currents  and  voltages  in  terns  of  the  state  variables  and  independent 
sourcea.  Thus,  the  rows  of  the  aatrices  A  and  B  are  cow posed  of  the 
coefficients  of  the  unique  linear  coobinations  of  the  state  variables 
and  the  independent  sources  which  are  equal  to  the  corresponding 
voltages  and  currents  on  the  left-hand  side. 

In  case  mutual  inductances  are  present  in  the  n-port  it  nay  still 
be  possible  to  choose  &  set  of  state  variables  as  specified  above  and 
write  Equations  (4-1);  however,  P  will  no  longer  be  a  diagonal  matrix. 
For  all  physically  realizable  n-ports  the  values  of  the  coefficients 
of  mutual  inductance  will  be  such  that  F  will  still  be  positive  definite 
and  symmetric. 

In  case  the  network  possesses  loops  which  contain  only  capacitors 
we  also  find  that  the  matrix  P  may  not  be  diagonal,  but  may  be  positive 
definite  and  symmetric  For  example,  for  the  network  of  Figure  4,1 


Figure  4 .  1 ..  A  network  with  a  capacitive  loop. 


3k 


we  may  choose  x(t)  <■  (v^(t),  v^(t)  )  ant*  wr^e 


<C1  *  v 

' 

x(t>  = 

0 

1 

O 

<C2  .  V 

0 

*1/5*2 

The  left-hand  side  of  this  equation  ia  equivalent  to  the  vector  (i.(t),ig(t) 
Clearly,  P  ia  positive  definite  and  symmetric.  A  similar  remark  can  be 
made  for  networks  having  cut  seta  which  contain  only  inductors. 

In  certain  of  the  cases  swntioned  above,  and  sometimes  when 
dependent  sources  are  present  in  the  network,  it  is  convenient  to 
choose  as  state  variables  linear  combinations  of  certain  branch  voltages 
and  currents.  In  any  event,  if  some  set  of  state  variables  may  be  chosen 
and  Equations  (4-1)  written  with  a  nonsingular  P  matrix  then,  upon  multiply¬ 
ing  both  aides  by  P*1,  we  obtain  an  equivalent  set  of  equations  in  the 
form  of  Equations  (4-1),  with  the  new  P  matrix  (the  identity  matrix) 
positive  definite  and  symmetric. 

In  addition  to  Equations  (4-1)  we  may  also  write 

w(t)  *  C  x(t)  ♦  D  u(t),  (4-2) 

where  w(t)  la  a  vector  whose  elements  are  the  remaining  port  variables 
(those  not  included  in  u(t)),  and  the  matrices  C  and  D  are  constructed 
in  such  a  manner  as  to  give  the  appropriate  linear  combinations  of 
state  variables  and  independent  sources  to  represent  these  port 
variables . 

If  we  consider  a  lumped  linear  multiport  netwsrk  containing  no 

independent  voltage  and  current  sources,  and  assume  thBt  at  each  port 

one  of  the  port  variables  (the  port  voltage  and  the  current  into  the 

port)  is  specified  independently,  then  it  may  be  possible  to  write 
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Equation#  (b-1)  and  (b-2)  us  described  above,  where  u(t)  Is  the 
vector  whose  components  are  the  independent  port  variables  and 
w(t)  is  the  vector  containing  the  remaining  port  variables.  Taking 
the  Lapiace  tranaTorm  of  Equations  (b-i)  and  (k-2)  we  easily  obtnin 

«{*)  -  [C(aP  -  A)'1  B  ♦  D)  U(e). 

At  fi  «  oo, 

W(s)  -  DU( a) • 

Thus,  if  we  consider  the  input  variables  to  the  multiport  to  be  +he 
components  of  the  vector  u  and  the  output  variables  to  be  the 
components  of  the  vector  v,  the  matrix  D  is  the  transmission  matrix 
lor  the  multiport  when  all  of  its  capacitors  are  short  circuited  and  all 
of  its  inductors  are  open  circuited.  If,  for  every  pair  of  diBtinct 
ports  there  exists  a  zero  of  transmission  at  s  ■  a,  then  D  is  a 
diagonal  matrix. 

(_■ .  A  Class  of  Nonlinear  Distr ibuted  Networks 

Many  nonlinear  distributed  networks  may  be  represented  as  in 
Figu.e  J*-2.  This  network  consists  of  three  main  parts:  One  part  is  a 
lumped  linear  multiport  which  is  connected  to  each  of  the  other  parts 
only  at  its  ports.  The  second  part  consists  of  the  collection  of 
lossless  transmission  lines.  This  part  is  divided  into  three  groups 
us  explained  below.  Each  end  of  each  line  is  connected  to  one  of  the 
ports  .!’  the  linear  multiport.  The  remaining  ports  of  the  linear  multi¬ 
part  are  connected  to  the  third  part  of  the  networx,  a  nonlinear 
;.r.i  1 1 1 ;  .  i  t. 
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Figure  4.2.  A  typical  nonlinear  distributed  network. 

The  nonlinear  Multiport  is  characterised  as  follows:  We  suppose 

that  there  are  nQ  ports  for  which  the  voltage  is  the  independent 

variable,  and  n.  ports  for  which  the  current  is  the  independent  variable. 
P 

and  let  v.  and  I_  denote  vectors  whose  coaponenta  are  the  values  of  these 
a  p 

port  variables.  Then,  the  remaining  port  variables  are  specified  by 
a  nonlinear  function  3: 

(£)-,<(£))-  <k-J) 

We  asauae  that: 

(Al)  It  is  possible  to  specify  vectors  u  and  v,  u  being  a  vector  whose 
components  consist  of  one  port  variable  from  each  port  of  the  lumped 
linear  multiport,  including  the  components  of  IQ  and  v^,  and  w  being  a 
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vector  whose  component*  consist  of  ths  reetaining  port  variables,  snd  s 
stats  vector  x  such  that  the  resulting  linear  network  may  be  characterised 
by  Equations  (m-i)  and  (^2),  with  P  a  positive  definite  syoawtric 
matrix,  and  such  that  there  exists  a  zero  of  trsnssission  at  t  •  *»  frost 
each  port  of  the  lumped  linear  multiport  to  any  other  port  at  which  a 
transmission  line  U  connected,  and  vice  versa. 

For  a  given  network,  in  order  that  the  lumped  linear  multiport  have 
the  required  transmission  zeros  es  specified  in  condition  (Al)  above, 
it  may  be  neceaaary  that  at  certain  porta  a  specific  port  variable  be 
assumed  Independent.  The  choice  of  independent  port  variables  divides 
the  collection  of  transmission  lines  into  three  groups:  One  group 
contains  all  of  the  lines  which  have  both  ends  connected  to  porta  for 
which  the  current  into  the  lumped  linear  multiport  is  chosen  as  the 
independent  port  variable.  Another  group  contains  all  of  the  lines  which 
have  one  end  connected  to  a  port  for  which  the  current  into  the  multi- 
port  is  the  independent  port  variable,  and  the  other  end  connected  to 
a  port  for  which  the  port  voltage  is  the  independent  variable.  The 
third  group  contains  all  of  the  remaining  transmission  lines,  each  of 
whose  ends  is  connected  to  a  port  at  which  the  port  voltage  is  chosen 
as  the  independent  port  variable. 


in  Equations  and  where  and  ij,  denote  and  nr  vectors 

Uy  *  n, }  whose  components  are  the  independent  port  variables  at  those 

ports  of  the  iump.ed  linear  multiport  to  which  a  transmission  line  of  the 

first  group  is  connected,  etc.  If  *  is  an  n*vector  then  p  and  A  are 

n  *  n  matrices,  B  is  an  n  x  (n  ♦  n_  ♦  . . .  ♦  n„)  matrix,  C  is  an 

a  p  v 

(n^  ♦  n^  ♦  ...  +  ng)  x  n  matrix,  and  D  is  an  (nQ  ♦  *  ...  ♦  n@)  x 

(n  ♦  n„  ♦  ...  +  nQ)  matrix.  We  let  B  and  C  be  partitioned  in  the 

Q  P  y 

following  manner; 

B.  tBr!Bn],  Bj-  [Ba!B9I,  Bu.  iB,.!^  |B,|  B^Bg], 


C“  e 

- >  Ur 


where  is  an  n  x  (n^  ♦  n  )  matrix,  Bg  is  an  n  x  (n^+  ...  +  n^)  matrix, 
B^is  an  n  xn^  matrix,  etc.  Finally,  we  let  be  partitioned 

by  columns  as 


V  lBr.iBrJ'"iBr  J'  B!>  *  <Bs iB6J  •  •  •  !B6 

12  n ..  1  n_ 


3.  etc.  , 


and  Set  C  ,  C&,  . 


c<:  partitioned  by  rows  as 
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We  also  assume  that  the  matrix  D  has  the  following  form: 


where  DT  is  an  (n,  +  n„)  x  (n  +  na)  matrix,  D„  is  an  n  X  matrix. 
J  upup  I  7  I 

etc.  Condition  (Al)  above  specifies  that  D^,  ...»  Dg  be  diagonal 
matrices. 


Vie  also  assume  that; 

(A2)  Each  of  the  diagonal  elements  of  D^,  D^,  . Eg  is  a  nonnegative 
real  number)  and  for  each  transmission  line,  at  least  one  of  the  two 
diagonal  elements  of  the  D  matrix  which  correspond  to  the  ports  to  which 
the  line  is  connected  to  the  lumped  linear  multiport,  is  a  positive  real 
number . 


From  Equations  (4-1),  (4-2),  and  (4-3),  we  obtain 


«  3(C  ^  +  D  T 


(4-t) 


We  assume  that: 

(A3)  It  is  possible 


solve  Equation  (4-4) 


lor  the  vector 


-t.  t 

V  as 


./pile it  function  of  C,  v.  in  some  neighborhood  of  CjX  =  0.  That  is, 
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we  as suae 


in  some  neighborhood  r  x  =  0,  G  =  (x  :  xe  E*1,  |jx||  <  H), 
that  there  exists  son  runetion  3*  such  that 

|  -  •«:£>.  (4-5) 

We  further  assume  that: 

(A4)  Equation  (4-3)  satisfies  a  Lipschitz  condition  in  G  ,  and  3*(0)  *  0. 

H 

If  conditions  (Al)  through  (A4)  above  are  satisfied  for  a  given 
nonlinear  distributed  network  then  it  is  a  member  of  the  class  of  networks 
for  which  our  stability  theory  applies.  It  is  felt  that  most  nonlinear 
distributed  networks  that  one  is  likely  to  encounter  will  satisfy  the 
above  four  conditions.  If,  however,  a  given  network  fails  to  satisfy 
one  or  more  of  these  conditions,  the  following  techniques  might  still 
be  used  to  render  it  amenable  to  the  application  of  our  theory. 

If  D  is  not  a  block  diagonal  matrix  with  diagonal  D^_,  D^,  . .  , 
submatrices  it  may  often  be  consistent  with  physical  reality  to 
consider  the  presence  of  small  "stray"  reactances  at  the  ports  of  the 
lumped  linear  multiport.  These  reactances  will  have  the  effect  of 
giving  the  necessary  zeros  of  transmission  at  s  =  ».  The  addition  of 
small  stray  reactances  at  those  ports  where  the  lumped  linear  and  nonlinear 
multiports  are  connected  will  always  allow  condition  (A3)  to  be  satisfied; 
for,  by  adding  enough  strays,  the  matrix  may  be  made  to  contain  all 
zeros,  and  hence  3*  s  3.  Finally,  if  the  function  3*  does  not  satisfy 
the  required  Lipschitz  condition  it  might  be  satisfactory  to  approximate 
3*  by  some  function  vhich  does--a  polynomial,  perhaps. 

In  our  theory  we  consider  only  undriven  nonlinear  distributed 

networks;  that  is,  the  networks  are  assumed  to  contain  no  independent 
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vo,  >ge  and  current  sources.  All  of  our  networks  have  an  equilibrium 
state  x(t)  *  0,  and  it  is  the  stability  of  this  equilibiium  which  we 
study.  If  a  network  contains  independent  sources  which  are  constant  for 
all  time  (bias  voltages,  for  example),  and  if  it  has  an  equilibriua 
state  other  than  5,  it  may  still  be  possible  to  study  the  stability  of 
this  equilibriua  by  first  finding  an  “equivalent’'  network  with  a 
corresponding  equilibrium  at  x  =  0.  For  example,  in  Figure  4.3  we 
have  shown  such  network^  this  network  has  three  equilibrium  points, 


Figure  4.3.  A  nonlinear  distributed  network  with  a  bias 
voltage. 

labeled  (a),  (b),  and  (c).  We  may  study  the  stability  of  any  of  these 
equilibria  by  considering  the  equivalent  network,  Figure  4.3a  with 
E  -  0  and  nonlinear  function  f  described  by  the  curve  of  Figure  4.3b 
with  the  origin  of  the  v-i  coordinates  shifted  to  the  appropriate  point 
^either  (a),  (b),  or  (c)). 

We  shall  now  derive  the  system  of  functional-differential  equations 
wnioh  describes  the  behavior  of  any  network  in  our  class 
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For  any  nonlinear  distributed  network  in  the  class  described  above 


we  have,  from  Equations  (4-2),  for  J  »  1,  2,  ry  =  n^. 


v„  (t)  -  C  x(t)  ♦  d  i  (t), 

h  rj  rj  ri 

and 

v.  (t)  =  C.  x(t)  +  df  i  (t). 

j 

Thus,  we  can  represent  our  network,  as  far  as  the  behavior  of  the 

j-th  transmission  line  connecting  the  r  and  £  ports  is  concerned,  by 

the  network  of  Figure  4.4a.  Hence,  since  0  <  d  <  °°  and/or  0  <  d.  <  », 

Tj 

we  have  (using  the  functional  equations  for  the  transmission  line  two-ports) 


(c) 

Figur*-  Equivalent  networks  onnected  to  typical  transmission 
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1  (t)  -  C  x(t)  *  2v  Ye*  C  x(t  -  2(lc»l)T  ) 

rJ  TlrJ  T.c0  rJ  Ti  rJ 


where  we  have  defined 


V(t)  =  c6.x(t)  +  VV10' 

<3  J  j 


kk 


and 


\(t)  ■  v(t)  *  %  v“ 


Thja,  we  can  represent,  our  network,  as  far  as  the  behavior  of  the  J-th 
transmission  line  connecting  the  6  and  T)  ports  ia  concerned  by  the  network 
of  Figure  4.4b.  Hence,  since  0  <  <  »  and/or  0  <  d^  <  •»,  we  have 

00 

l5  <*>  a  S  C6  *(t)  *  2%  Z  *1  Cft  *(t  '  2<*+l>V> 


3  J 


J  R«0  ^ 


av  I  o»,V(‘  - 


k-0 


j  'J 


w 

%  .<*>  -  -\.®,  «<•■>  *  I  VV*<‘  -  *C»U«,  > 


\  ') 


k»0 


b  "a 


*Y  I  -  t2**1*,  > 


^  ^  ’j 
J  k-0  J  J 


where  we  have  defined 


V  ‘  a.  *  z  ' 

J  ’j  % 

J 

z  r 

°6ni  '3 

\  '  (d  +  z0  )2' 
J  6lb 


n  d  +  Y  * 

J  "j  °">j 

-Y  l’ 

°6<l,  6J 

M  =  - J - J"  , 

j  (d  *  Y  )2 

nj  Bn, 


^  ^  ’  i\  *  \\  *  v 
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and 


Thus,  if  we  lee  a.  fr^,  ....  r„  ,  6,..  B.,, 
end  define  S'  .  !'(i)  on  B  a»  follovd: 


t'  •  for  |  .  rv,  0  -  1,  2,  ....  nr=  n{), 
!'  =  nj  for  5-6.  0-1,  2,  ....  ^  .  n^}, 

{-  -  for  {  .  Cj  (j  .  1,  2,  ....  „f  .  „e), 

l'  ‘ri f°r !  *  5j  <J  -  •••-  "t  -  v. 

!'  "  »j  for  i  =  1j  (J  -  1,  2,  ....  n^  . 

•'  ‘  £j  r“  ♦  ■  ®J  (J  ■  1.  2 . n6  .  nj, 

we  have 

W-IV 


2  £  ptMtBiV(t  ‘  2{k+1)V 


Jt»0 


2  ZPt  VJBtC{'x(t  *  )]. 


k-0 

Substituting  the  above  equation  and  Equation  (4-^)  into  Equation  (4-1), 
we  obtain 


F*(t)  -  Ax(t)  ♦  Bj  fl*(Cj*(t))  -  Y, 

{£U 

00 

*  I  »h;<‘  - 

(cQ  k-0 
00 

*  2Z  Z  -  <2k+1)T{)> 

(€U  k-0 

where  we  have  defined,  for  Jcfl, 


(4-6) 


Q  •  U  C  ,  M  :  |i  9  C  i  N  •  v  B  C  .  . 

t  tit1  i  i  t  «'  t  tit 

Multiplying  each  aide  of  Equation  (4-6)  by  P  1  (which  exists  since  P 
is  positive  definite,  by  condition  (Al))  yields  a  functional-differential 
equation  of  the  desired  form  which  describes  the  behavior  of  the  network. 


4 .  The  Llpschitz  Condition 

It  is  easily  shown  that  the  right-hand  side  of  the  functional-differential 
equation  which  describes  the  behavior  of  any  network  in  our  class  satisfies 
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I A  |A 


a  Lipschitz  condition  on  C^:  If  we  let  xfc  be  a  point  in  C^,  then  by 
Leona  1  of  Section  the  sapping  f  froo  Cu  io  I  defined  by  7 (x  }  m 
x,(0)  =  x(t),  satisfies  a  Lipschitz  condition.  Thus,  so  does  the  napping 

v 

7  defined  by  f,(x. )  ■  P  ^Axft),  since  if  there  exists  L  such  that  }|${0)-f{0)|| 

Lp(i,t)  for  every  q,  *  •  C^,  then  -  P*  1A*<0)(|  <  f|P  1A)| *  jj^( 0 ) 0 > ji 

UlP~1A|!-t/)p(5i*)-  Alao>  the  sapping  7(xt)  s  tf»(cj(t))  satisfies  a 

Lipschitz  condition,  by  condition  (A4),  and  hence  so  does  the  sapping 

f  (xt)  a  p"^ ^*(Cj x(t)) .  Finally  for  each  l,  the  sapping  7^  of  CH  into 

E  defined  by 
n 


7  fij  =  -P’^xtt) 


2P’Hl  x(t  -  (2K+1)t4> 


has  the  form  of  Equation  (3-12),  which  has  been  shown  to  satisfy  a 
Lipschitz  condition.  Thus,  a  simple  application  of  the  triangle  inequality 


shows  that  the  mapping 


T  +  £ 


vhich  is  the  right-hand 


side  of  our  functional-differential  equation,  satisfies  a  Lipschitz  condition. 


Chapter  V 

STABILITY  (V  NOXUMlAft  DISTRIBUTE)  MWWOWffi 

It  ia  obvioua  that  *  •  0  la  an  equilibrium  solution  of  the  functional - 
differential  equation  which  describes  the  behavior  of  any  network  in  the 
claaa  of  nonlinear  die tribute*  networks  define*  in  the  previous  chapter. 

In  this  chapter  we  state  and  prove  several  theorem*  concerning  the 
stability  and  instability  of  thia  equilibrium  solution.  We  use  the  same 
notation  as  In  Chapter  IV  and  consider  our  distributed  networks  to  be 
characterised  by  Equation  (U-6),  with  €  C^. 

1.  A  Lyapunov  Functional 

Before  defining  our  Iyapunov  functional,  V,  we  shall  prove  several 
useful  Leonas. 

Lemma  l.  If  A  is  an  n  x  a  matrix  and  if  B  -  AAt,  then  ||B||  -  O  if 
and  only  if  A  ■  0. 

Proof:  ("If")  Proof  of  this  part  is  trivial.  ("Only  if")  Since  B  is 
symmetric  there  exists  a  nonaingular  matrix  P  such  that  P  -1BP  ■  A,  where 
A. is  a  diagonal  matrix  whose  diagonal  elements  are  the  eigenvalues  of  B. 
Since  B  is  positive  semidefinite  all  elements  of  the  main  diagonal  are 
nonnegative;  but  since  (|B||  ■  0,  all  elements  of  the  main  diagonal  are, 
in  fact,  zero.  Thus,  A  ■  0  and  hence,  B  *  PAP  1  *  0.  If  Bj.  represents 
the  element  of  B  in  the  i-th  row  and  j-th  column,  and  similarly  for  A, 
m 

"  Y.  Ai*v 

k*l 

In  particular,  for  i  «  1,  . . . ,  n,  ■> 

AiR  -  0.  Thus,  A  *  0.  Q.I.D. 
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Aik  and  hence,  for  k 


1. 


m, 


k»  1 


Leaaa  2 ■  If  A,  B,  and  p  art?  any  real  numbers  such  that  A  >  0  and 


jp|  <  1,  and  x  and  y  are  arbitrary  n-vectors,  then  the  following  inequality 
holds; 

PS  pk<x,y)  -Ajpj*<y,y>  <  |  p|K  <  x,  x  )  ,  for  k  =  0,  1,2,  ... 

Proof;  By  the  Schwarz.  inequality, 

|<*,y>|  <  lixIMlylU  *  *,  y  ‘  En. 

But  for  any  x,  y,  (j|xj|  -  ||y||)c  >0  and  hence, 


rf  *  t|y||2  >  2|WMW|. 

Therefore, 

2\  <x,y  >  |  <  !|x|f  '  !!yl!2 

and  hence, 

L  2<x,y)  <  <x,x>  ♦  <y,y> 

-* 

+  ;'lf:!k<*.y>  i  Mk<x»>0  ♦ 

|f>lk<y,y; 

1  •'|^Jk<x,y/  -  t>'!k<y,y;  < 

M*  <  x,x) 

- 

2^  >.  x.y  -  |  p | k  <.  y , y  ,1  <  |p| 

ik  / - . 

\x,;</ 

.’i~k  <  Bx.Ay  -  |  ?-lk  Ay,  Ay  <  |  uj  K  <  Bx,  Bx  ✓ 

.. . .  k -  „  2  .  .  k - v 

.Are  x ,  y  ,  -  A  |u|  \y,y> 

1  !^|\x.x, 

B*  *  <  x.y  ,  *  Ajt  jk  <y,y  ,  v 

'/  ,  ,  k - 

^  !  *  |  ^  A  >  X  x  ■ 

Lem  a  ' 

1*  A,  :* *  'j.'.i  x  m*  'ii,y  r*  il 

•  r ..  .'  *  u  ‘ 

s  ana  x  a.'.J  y  if-  art. trary  ».«*.•••  •* 
i.-  ■ 
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NX 


AjoAy.jr>  »  *  ~~ 

ik 

for  K  =■  0,  1 ,  2,  ... 

Proo  f :  With  only  a  few  sign  changes,  the  proof  is  identical  to  that 
of  T.-emm*  ?■ 

We  now  define  a  Lyapunov  functional  V  on  C*  *  U  C  .  For  each  JcU 

Q<X<u>T 

let  a^  and  denote  fixed  nonnegative  real  numbers,  the  values  of  which 
will  be  chosen  later.  For  }j  >  0,  define  the  functional  V  for  every  <p 
in  C„(  (—.0],  E11  )  by 


V(?)  =  \  *\o)  P  9(0)  ♦ 


-2ki. 


iZ!  !p|lk  I  £pt(p)‘p(o)dt 


{cil  '  k=0  (7—  2(k+l)t 


i 


-(2k+lk 


+b, 


r 

|  /V'C  0)9(0)  do  +  Y-Ip|  1  f 

'  ~  o=-(2(k*l)H)T 


^(c)  9(0)  do 


o=- 1 


k=0 


to-u 


Along  a  trajectory  of  our  system'  we  have 

(«  t-2ki^ 

j  ^(o)  *(c}  a° 

k=0  . 


c=t-2(k+l)i 


t- (2k+ 1)  t 


* 


j  xt(0)  x(o)  do  +£j0^k+1  /  **(0)  *(<>) 


0=t-t 


l 


k=0  o=t-(2(k<-i)+l)r 


i 


The  derivative  of  this  functional  along  a  trajectory  of  the  system  has 
the  value 


37  V(xt)  *  ^(t)  P  x(t) 


■lid  k 

$tfl  '  k=0 


I  f  xt(t-2ki^)  x(t-2kt^)  - 
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-  *t(t  -  t^)  x(t'T^)  ♦  -  (2k*l)i^  )x(t-(2k*l)‘r^) 

K=0 


-x^t  -  (2(k+l)+l  )t  )  x(t  -  (?(k+ 1 ) 


1)^)1 


where  P  x(t)  may  be  replaced  by  the  right  hand  side  of  Equation  (i*-6). 
Before  making  this  substitution,  however,  we  rewrite  this  equation  as 


d_ 

dt 


V(xt)  =  xL(t)  Px(t)  + 

*fC 


I  *t  [  Z  'v  “  ;v 


t  eU  k=  1 


2k -r^  )  x(t  -  2k -r^  ) 


“  1 

-  £  |p^|k  Xt(t  -  2(k+l)^)  x(t  -  2(k+l)i^)| 

k=0  J 

X,  b  x*(t)  x(t)  -  ]T  xl(t  -  i  )  x(t  -  t*  ) 


tea 


tea 


]C  bt  X!let^K+1  yt(%  *  )x(t  -  (2k+l)T^) 


t  «.a  k~o 


00 

-  Xj0!1*  xb(t  -  (2k«-l)Tft)5(t-(2k* 1)^)1 
k=l  J 

x^iiP  x(t)  •*  ^  a^xb(t)x(t)  *  ^  x^(t)  x(t) 

5  <a  t 

00 

^  >t(t-?(k+l)TJ ;  )x(t-2(k+l)Tj  ) 


f  fii 


k=0 


^  (  Ip^I  -1)^  *0|l  k  Xt(t-(2kt  :  ytj  >x(  1- t-  ’K-e  1 )  )  . 


f  *a 


k*0 
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a.  if* 


V(xt)  -  ^(t)  Ax(t)  ♦  t)  Bj3*(Cjx(t) )  -  xl(t)  Zs  x(t) 

ir  a 

*  ^(tj  *  ^(t)£  b{Inx(t) 

4  cu  4 cu 

+  Z  Z  [2p4  *t(t)  M4 

4efl  k=0  L 

-a^(l-|p4l  )  lP|l  K  x*(t  -  2(k>l)i5)  x(t  -  2(R+1)t6', 

♦  £  £][2P|  ^(t)  H&  x(t  -  (2k*l)t  ) 
tea  k*=cr 

-bt(l  -|ptl  )  lp6l  k  xHt  -  (2k+l)t4 )x(t  -  (2fc+l)t{)  , 


V(xt)  *  xt(t)Ax(t)  «■  x^t)  BI3*(CIx(t))  -  Q|x(t) 

4*2 

+  x^t)  £  a^Inx{t)  ♦  *t,£  b^n  x(b) 
tea  | eC 

00 

+  £  £  *(■>»  *(*  •  2(k+l)i| )  > 

4  eU  k*0  L 

-a  (l-|p?|  )|0|lK<x(t  -  2(KU)t4),  x(t  -  ?(k+  1)t|  )  > 

+  T.  Z!  [2p4 '  Mt  *(  1 ’  *(  *  “  (2k+ 1  h4  }  ^ 

4rfi  k^OL 

-bt(l-|pel  )|Pf|  \“(-'  '  (?kM)Tt),  x(t  -  (2K*1)t4)> 
5** 


In  cam  M.  0  !Vr  -r»o  »f  £,  let  us-  -  O;  and  In  case  R  •  0 

for  some  Ipt  u.i  pi  -h  a  0.  Otherwise,  we  shall  require  a^  >  0 

and  b,  >  O.  That  this  is  an  appropriate  choice  of  a  and  b  for  these 
3  5  £ 

cases  ran  be  by  comparing  Equation  (4-ti)  vith  ,ur  Lyapunov  functional 

V  in  Equation  (}-]).  If  H,  *  0  then  the  terms  involving  x(t  -  2(k+l)t^) 

ar*  not  present  in  Equation  (4-6).  Similarly,  if  N  =  0,  the  terms 

> 

involving  x(t  -  (2K+l)if)  are  not  present  in  Equation  (4-6).  Thus,  for 
such  values  of  |  we  have  no  need  for  the  terms 


Z>/ jv 


(o)  q(o)  do, 


<?=-2(k+  1)t, 


0  •  -(2k+l)^ 

j  4t(o)  4(o)  do  ♦  IPjl^1  J  5t(o)^(o)  do 

ofc-i.  k*0  o=.-(2(k+l)+l  )t- 


in  our  Lyapunov  functional.  Vie  also  note,  by  Lemma  1,  that  *  0  if  and 
only  if  =  0;  and  -  0  if  and  only  if  ||N^N^|j  *  0.  To  avoid  awkward 

notation  we  shull  use  the  following  conventions:  the  symbols 


denote  summat.cns  over  only  those  {tii  for  wnich  /  0  and  b^  f  0, 


1 espectively .  Thus, 


J[V(xt)  =  xt(*  1  fiJl)  1  v:(t)  BIS5*(Clx(t,) 


s,x(,/  •  ;t(t,E  a?V''y  4  b».v(t; 
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,  ® 

n  i.  _  .r*  * 


{IU  k=0 


“as(i*  OU  -  2(k*i}t^},  *(t  -  i?(k+ 1  }T| )  > 


£  11  [2°t  *0<  -  (Pk+l)tj )  > 


tea  k*0 


-b|(l  -IpjI  )|p6|k<x(t  -  (2k+l)T{),  x(t  -  (2Rn)ie)> 


But  then,  using  Lemma  2,  we  have 


ft  V(*t)  <  ^(t)  Ax(t)  «  *t(t)BI»»(C^(t)) 


-*t(t)  £<^x(t)  *  ^(t)£  afcInx(t)  ♦  *t)£  bftXnx(+.) 

4*  a  810  US 

/  «» 

*Z  Z  ..a1-  Mt  ;(t)> 

*=0  5 

✓  /  00 

*  Z  I  ETrrqrn  *£<■''<>  ■ 


it  a  k=o 


£  Ip, |K  -  i/U  -  |pj|), 


wo  have 


~  V(xt)  <•  xl(i)  Ax(t)  *  xl(t.)  Bt3*(CEx(t^ J 


-x  C) 


Z  s*(t) 


“t  o 

*  (O  /  (a  I  ♦  - L _  y  |(t  ]  ^ 

£  *“  IP,!)'  W1*'1’ 


•  «V)  l 


L  *  -77-7-73  Yi1  <(l>' 
e*  u  b{  ( i  -  I  ^  I )  ?  ? 


01  v(*t}  -  0 


for  all  7t  C  CH  If 


(•1^- •  ,£-  [>(1] 


U4(J  *  1^1) 


~*s*Stj+F  £b. i  ♦ — 

5  t  n  .  ,  2  N*M»  J 1  x> 

etu  v1  -  |*'4| )  5  5 


for  .11  ;  .  *»,  fin  <  «  i.  We  nottt  UMn  ttmi 


SJn  f  ”  j  772  +  y  [b  I  +  —  1  1 _ _  M  - 

*to-i»si)^  (t;  »"  ita-  |p^s  V V* 


==  <r  • 


V>|u( 


—2  IIHjNJii]]  •  Hxir ; 


r  VKJ  *-  0 


ir  all  x,  C  CH  if 
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i 


(“u 


y 


th.  1 


~  sw 


V  lb  • - 1 - ;DN.».t||], 

r«  v-k^  !! 


for  all  x  C.  Gu,  x  /  0.  In  order  to  obtain  conditions  for  ensuring  that 
H 

^  V(x^)  <  o  that  are  in  general  us  weak  aa  possible  for  this  method,  wc 
shall  now  choose  the  constant,;;  /  0  and  b^  >  0  such  that 


u.  * 


1  atU-!ptir 


and 


[b.  + 


Ilyvll  ) 


1  bf(l-in,l)"  1 

-*  » 

are  minimized.  The  funs.' lion  u(.i)  -  s  +  k/s,  k  >  0,  has  a  minimum  for 
c  >  0  at  s  =  the T" to r e ,  it  is  <  lear  that  we  should  choose 


\  =  ~u 


Vifor 


it-.-i 


V  "*>  *h  " 
i  -  -> — * — 5 - 

«  o  1  u 


for  t  t  a,  Uj,  4  0; 


V  r  t  t  a,  b  40. 


, :  ’  h  this  i.-hc  :  •  i  and  t  we  may  writ'.' 

*  r 


^va, xivtiu-  .y  ji'Vi11,', 

«’•  1  ~  L*  u  *  ip.  i  .* 


St  M  Si  M 


fl  +  - : - Kll) 

,pt^  '  "  5  * 


♦  )  — 5"i  (  I  ♦  - K - N.  N 

1 —  U  "  |=?|  ■  '  n  (tv  vln  4  ] 


4t« 


*Vi» 


Vi  *('.) 


♦  xV)  B  B*(C  x(t)). 


Defining 


A  = 


"  v  O 

O  V 


(5*2) 


we  note  that 


I5!  '  EVi5!  'Z* t\C(  *  lBr  IVi-V  1 

fi-.  SCU  S6w  1  ‘  nB 


’  [Br  'Rr  !•■•!*«  !  A 

i  -•  na 


ne 


Vcr 

i  i 

K.a 

V>  r0 


\  c* 

nB  n© 


B  A  c- 
C  u 
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w*>  now  dpfir\»‘,  *, 


iti\  :i  K 


u. 


tr  lift  ii  i  <> 

* 


if  EJb^li  «  v 


K„  »  (K  >E  »K.  1, 

U  Ij  I  v  *  v 


r  tii'n 


^-i '  yil^^  li 


$C  u 


rr?IVT  \ 


Z;  - 1^1)  \  ^  itceir?  iib^ir-  1  ? 


1m.  |  •  r.  ii  •  hb^ 


<r'  ‘"t1 

a  l . a  *•  i;;r 


A-  re 

a  i 


t  t.ii 


f  uii 

add  «imi  lar  ly, 


V-  ^HiV1  i:v  _.t 
b  . \ n  \ 


r- « *  /tN.  N»  I!  f  . 

£  l  ^7  V. 


>,  MV  INIb  II 

a  ,,. i-  V  V  V 


ri  ",  l  H 

t  u  * 


it'  Wt*  dfrf i fl*.‘ 


dl 


U^Nlc|ll  ♦  I vt I  ‘ IfCj,  il ) ■  I|b^|| 
u  -  Ip. D 


for  f.  <  U, 


o 

O  " 


(5-3) 


we  have 


ft  V<*t>  ^  *'(t)  [A  -  Bc  A  CD  ♦  En  tE*  ♦  (tr  f)In]x(t) 


♦  xL(t)  BI3*(CJ.x(t)) 


We  may  also  define  a  functional  V  on  C*  by 


I  “  r  "21tT| 

v(?)  =  ^  ^(O)p  ?(0)  |p^l  Ky^t(°)  ?(°)  do 

\  k=0  os-crCk+l)^ 


r-  y  00  M 

+  b^l  /  do  +  I °4 1 K+  J  ^(o)  do|  )  ' 

0=-t5  *=G  o— (2(k*l)+l)i6  *(>5) 

Equation  (5*5)  is  identical  to  Equation  (5-1)  except  thrt  the  first  plus 
sign  is  changed  to  a  minus  sign.  Proceeding  in  exactly  the  same  manner 
as  above,  with  the  use  of  Lemnu  5  instead  of  Lemma  L‘,  one  easily  obtains: 

^  V(xL)  >  7t(t)[A  -  Bn  A  Ca  -  Ku  fEj-  ( .r  'OlJ  x(u 

+  x*(t)  J(t) ) .  [•>-(■) 


.-(2k+l)i. 


ienotes  the  tract'  of  fht-  matrix  i. 


One  may  easily  show  that  for  any  H  >  0  the  functions  V  of  Equation*  (5*1) 


and  (5-5)  are  continuous  and  bounded  on  Cyj  this  ia  done  in  Appendix  C. 

It  is  also  obvious  that  7  of  Equation  (5*l)  satisfies  V(^)  >  0  for 
at)  «f  C*. 

We  now  state  and  prove  the  basic  stability  and  instability  theorems 
for  nonlinear  distributed  networks. 

2.  Stability  and  Instability  Theorems 

The  following  theorems  apply  to  any  network  in  the  class  of  nonlinear 
distributed  networks  defined  in  the  previous  chapter.  That  is,  we  assume 
that  conditions  (Al)  through  (A4)  of  Section  4-2  are  satisfied  by  every 
network  for  which  these  theorems  are  to  be  used. 

Theorem  1.  If 


K  (2V*p)^2>  where  kp  is  the  smal’sst  eigenvalue  of  the  matrix  P, 
and  if  xHa  -  BQ  A  CQ  +  En  V  Eg  +  (tr  lf)l  Jx  ♦  <  0  for 

all  x  in  a  neighborhood  of  the  origin  GR  «  (x  :  x  €  E1*,  ||x||  <  K), 
then  the  solution  x  *  O  of  Equation  (4-6)  is  stable.  Furthermore,  if 
xt[A  -  B  A  Cu  +  Eflf  E*  +  (tr  *)!  lx  ♦  xS-  d*(Cp)  <  0  for  all 
x  t  G^,  x  ^  0,  then  the  solution  x  =  0  is  asymptotically  stable  and 
every  solution  of  Equation  (4-6)  with  initial  condition  xq  in  Cy  approaches 
zero  as  t—*-». 

Proof.  Let  the  continuous  functional  V  be  defined  on  C *  by  Equation  (5-1) • 
Clearly,  V(0)  -  0.  For  the  riven  /,  let  denote  that  region  of  C* 


d3 


where  v(£)  <  i.  it  ut-!“  V(i)  "  < lli(o)f  <  ^(O)?  V(i)  <  i 

|!<p(0)|i  <  (Si/^p)^?  I!.  Let  u(s)  be  defined  on  (0,K)  by  u(s)  =  ijjApS2. 
Clearly,  u(s)  i*  continuous  and  increasing  for  0  <  s  c  K,  and  u(0)  =  0. 

From  Appendix  C  we  find  that  q  £  CH«^V(cp)  <  f,  and  hence  CH  c  U^.  Clearly, 
u(ll<p(0)||)  -  |y<p(0)f  <  V(<p),  for  all  q>  G  U,,  as  observed  above.  According 
to  Equation  (5-4),  the  condition  that  x^lA  *  Bjj  A  Cg  +  ¥  E^j  +  (tr 

*  xtBitJ*(CIx)  <  0  for  all  x  C  GK  implies  that  <  0  for  all  q>  £  U^. 

Hence,  by  Theorem  2  of  Section  2-2  the  solution  x  =  0  of  Equation  (4-6) 
is  stable. 

The  condition  that  x^A  -  A  Cu  +  EQ  f  +  (tr  ¥)In]x  ♦  xSj  3*(CIx) 

<  0  for  all  x  C  G^,  x  ^  0,  implies  that  M,  the  largest  invariant  set  in 

R  (the  set  of  all  points  in  where  =  0),  contains  only  the 

-  -  2 

point  <p  =  0.  Thus  by  Theorem  1  of  Section  2-2,  every  solution  of  Equation 

(4-6)  with  initial  condition  x^  in  U.  (in  particular,  all  x  £  C„)  approaches 

o  l  on 

zero  as  t-*«e.  The  solution  x  =  0  is  therefore  asymptotically  stable.  Q.E.D. 

Theorem  2.  If  x^A  -  B u  A  C Q  -  E3  ¥  e\  -  (tr  ¥)I  ]x  ♦  xSt  3*(CIx )  >  0 

for  all  x  /  0  in  some  neighborhood  of  the  origin,  G^.  =  (x  :  x  G  En, 

IJx||  <  ^  H),  then  the  solution  x  -  0  of  Equation  (4-6)  is  unstable. 

Proof.  Let  the  bounded  continuous  functional  V  be  defined  on  Cj,  by 
Equation  (5-5)*  Let  Uy  denote  that  region  of  Cy  where  <p^(0)  P cp(0)  > 

^See  pp.  110-111  of  reference  [37]- 

This  is  seen  as  follows:  Clearly  {0]  is  an  invariant  set  in  R.  Now 
if  5  G  R  then  q(0)  =  0  (since  q  t  ,  q(0)  0  =*  V^_^j((p)  <  0).  Suppose 
q  c  M;  then,  according  to  the  definition  of  an  invariant  set,  3  a  function 
x,  defined  on  (-«,«),  with  v.  M  v  t  in  (-»,«)  and  x^  =  q>.  But  then, 

for  all  t  ^  0i  x^  £  M  c  R  —  x^(O')  ■  Q.  That  is,  x(t  )  -  0  for  -<»  <  t  <  0. 

Hence,  <p  *  0. 
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/  *  ‘PkT»  r  0 

y [ay  yKr/ ^t(-)  ♦  bi[y  d° 


4<.u'  kaO  0=-£?(k*  l)  t 


l 


<^-t. 


*  -  (2K+  1 }  -i^ 

y  ipjj k*iy?t(o)  i(o)  do 

k=0  <*-(l>(k+l)U)i. 


that  is,  denotes  that  region  of  where  V(q>)  >  0.  Now,  since  (z  : 

•/.  e\  2  >0)  is  open  in  E1,  and  since  V  is  a  continuous  mapping  of  C^ 

into  e\  U  is  open  in  C^..  Thus,  3  an  open  set  U  in  C  such  that 

V^cp)  s  0  on  U  =  11  Pi  Cy. 1  Clearly,  that  part  of  the  boundary  of  U, 

c»U,  which  is  in  C  consists  of  the  collection  of  all  those  points  <p  in 

C  for  which  V(q>)  -  0.  It  is  also  clear  that  5  belongs  to  the  closure 

of  U  =  l?  D  C^,  since  0  C  C^,  and  V(0)  c  0^0  fc  diy.  Let  us  define  the 

function  u(s)  mapping  the  interval  [0,H)  into  E  by  u(s)  a  -\pS' ,  where 

kp  is  the  largest  eigenvalue  of  the  matrix  P.  Clearly,  u(s)  is  continuous 

and  increasing  for  0  <  s  <  H  and,  u(0)  «  0.  Furthermore,  u(||<p(0)||)  °- 

^p||q(0)|l2  >  |  ^(0)?  cp(0).  Hence,  V({ji)  <  u(||i(0)||)  on  Uy*  U  fl  C  - 

According  to  Equation  (5-6),  the  condition  that  xt[A  -  A  C^  -  f  • 

1  tr  i )  Jx  +  tl*(C,x)  >  0  <r  x  <■_  G^,  x  ^  0,  implies  that  >  0 

on  the  closure  of  U  -  I)  d  C  ,  and  that  the  set  of  o  in  the  closure  of  U 
r  v*  ^  r 

su:h  that  -  0  contains  no  invariant  set  of  Equation  (**-6)  except 

p  =  0.  Thus,  by  Theorem  ?  of  Section  2-2,  the  solution  x  *  0  of  Equation 
(s->.  )  is  unstable.  Q.  E.  P. 

Several  corollarl'  s  to  *:.*.•  above  theorems  may  also  be  stated.  One 
M-'lug  a  trivial  prod'  is: 

■  i 1  *  existence  of  the  open  •  U  fo i lows  from  Theorem  it., 

:  •  rt-nce  f.?2j. 
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Corollary  1  (Complete  Stability)-  II'  conditions  (A})  and  ( Aa )  of  Section 
1*-;?  hold  for  all  H  >  0  und,  if  xl[A  •  A  Cu  +  Efi  f  (tr  * 

♦  (C-pc)  <  0  for  all  x  ^  En,  x  /  0,  then  every  solution  of  Equation 

x  -l. 

\^-6)  with  bounded  initial  condition  ^  £  C  approaches  zero  as  t — 

and  the  system  is  asymptotical  1>t  stable  in  the  large  (completely  stable). 

Proof.  Let  x  be  some  hounded  initial  condition  in  C.  Then.  3  H  >  0 
-  o 

such  that  ||xo(t,)||  <  H  for  all  t,  £  (“«“,0]}  that  is  xq  t.  C^.  We  may  use 
this  Value  of  H  in  Theorem  1.  Q.E.l). 

Before  stating  the  remaining  corollaries  we  define  the  concept  of 
a  critical  point  of  a  mapping  f  from  E*1  to  e\.  and  state  a  well-known 
theorem  and  two  lemmas. 

If  f  is  a  differentiable  mapping  from  E°  to  El  and  if  for  xq  E°, 
of/dx^  (xQ)  =  0  for  i  *  l,  n,  then  xQ  is  said  to  be  a  critical  point 

of  f  A  theorem  which  is  available  in  many  references  [}8  p.  62,  59 
p.  6l]  is  the  following: 

Theorem.  If  f  Is  a  twice  continuously  differentiable  function  mapping 

a  neighborhood  N  c  En  of  x  into  E  ,  and  if  x  £  N  Is  a  critical  point 

o  o 

of  f  then,  if  the  quadratic  form 

n  n 

y  y  ••••••  -,v:  -(x )  w 

/  .  /  .  dx , Ox .  o  i  j 

l-l  ,1-1  !  •’ 

is  positive  definite,  !'  has  a  strict  relative  minimum  at  x^j  that  is, 
•tiff  "Xiats  an  open  set  <i  N  containing  x  such  that  c  G,  x  /  <=» 

■'Gi  >  m'x  ). 

u 

ire-  hypotheses  of  Theorems  i  and  2  of  this  section  involve  conditions 
’ tri  -'un  -tions  f  of  the  form  V-x)  ~  x*Qx  ♦  x^B  B*(C,  *) ,  where  Q  is  an 
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n  x  n  matrix,  be  nonnegntive  or  nonpositive  in  some  open  set  Gjj  contain¬ 
ing  x  =  0.  We  may  use  the.  above  theorem  to  give  conditions  which  may 
be  easier  to  verify.  We  first  prove  two  simple  lemmas: 

Lemma  1.  Let  f(x)  s  x'^x  ♦  3«( C^x).  If 


dx  ox 
*  J 


denotes  the  matrix  whose  i-jth  element  is 


<•»'"  Isfsj  (°>j  *  « *  V'<5)cr  ‘  «*  • 

where  3*^(0)  denotes  the  Jacobian  matrix  of  the  mapping  3*,  evaluated 


^  dXj 


Proof. 

df  dx*  .  — t  .  dx  ,  dx*'  „  — 1„  dx 

r —  -  5 —  Qx  +  x  Q  r —  +  r —  B,.3#(C_x)  +  x  B-3*  (C_x)C.,  r — 

ox .  dx,  v  %  dx  dx,  I  I  j  '  I  '  1  dx, 

4-  J  J  J  J 


v-t  s“  ■  -t  ■  “  v-t 

dx  ..  „  dx  .  dx  „  dx  ,  dx  „ 


dx  .  dx'' 


..a.f  -  (0)  =  i^Q2i-e2L.Q22-+S2-P  5*^(0}  c  -  +  ^  B  3*'fo)  C 

ox.  ox.  '  dx  M  dx,  ox,  *  dx ,  ox,  ±u  '  '  i  dx.  dx,  l"  '  1  ox. 

*  J  J  1  *  J  J  11  J 


■  «ji  *  «ij  *  te,3-'Co)c,JJ1  . 

■  <S  '  r-t.i  *  *  t*,**^)  . 


Hence, 


TTTT  !5) 

J 


-  Q  +  *  Q*  ♦  [BI3»/(;)  Cjl\  Q.E.D. 


'<*  ‘  A‘  *D  A  co  i  Ku  f  Eu  -  * r  ^ T„ ’ 
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Lemma  2.  The  matrix  M  ♦  is  positive  definite  if  am'  only  ir  M  io 
positive  definite. 

Proof.  ("If")  Proof  of  this  part  is  trivial. 

("Only  if")  Let  x^H  -*  M^x  >  0  .  Then,  xStx  ♦  xV'x  >0  -»  kxl  Hx  >  0 
x'Mx  >  0.  Q.E.L. 


Leittnas  1  arid  2  prove  that  if  f(x)  *»  S*(G^  x )  then  the 


matrix 


^  (0) 


Sx,  dx . 
i  0 


is  positive  definite  if  and  only  if  the  matrix  Q  +  v5'»/(0)Cj  Is  positive 

drfinite.  We  now  state  corollaries  of  Theorems  1  and  2: 


Corollary  2  (Asymptotic  Stability).  If,  for  all  x  in  some  open  set 
containing  the  origin,  G^  =  fx  ;  *C  En,  ||x|j  <  K ),  the  function  xSjB^CCjX) 
mapping  G^  into  E1  has  continuous  second  partial  derivatives,  and  if  the 
matrix  -A  ♦  B^A  ¥  E^  -  (tr  ’J' )  *  BTB**(0)Cj,  where  B^O)  denotes 

the  Jacobian  matrix  of  the  mapping  B*  evaluated  at  the  origin,  is  positive 
definite,  then  the  solution  x  ■  0  of  Equation  (i*-6)  is  asymptotically 


:  table. 


Proof.  Letting  f(x)  *  -xt[A  -  B„  A  C,  +  E.t  e!  +  (tr  11)1  ]x  -  x1  B.3»(C^t), 
-  U  a  0  □  n  11 

we  see  that  f(0)  *  0.  The  hypotheses  of  this  corollary  imply  that  the 


quadratic  form 

n  n 

E  E 

i--i  j*i 


Of  ,T,.  , 

- 1 -  (0)A  a. 

oxi  ix,  i  j 


is  positive  definite.  Clearly, 
fee  i.  »  1,  .  .  .,  n,  and  thus,  x 
t  h-  aoove  theorem,  there  exists 


(see  first  line  of  proof  of  Lemma  l)  of/ox^  - 
0  is  a  critical  point  of  f.  Hence,  by 


an  open  set  G^,  c  G^,  containing  x  -  0, 
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that  x  £  G„  , ,  x  /  ^  '  <(<)  >  f>«=*-f(x)  *.  0.  Thus,  by  Theorem  1, 

ft 

*h<  solution  x  -  0  oT  Equal  i-n  is  asympt.nti  cany  .liable.  Q.E.D. 


Corollary  3  ( Instab il i ty) .  If,  for  all  x  in  some  open  set  containing 

the  origin,  G  =  (x  :  x  <\  En,  ||x(|  <  y  1,  the  function  x*8  8*(C  x)  mapping 
T  l  x 

tl  into  E*  has  <-u»»Un««u4e.  second  partial  derivatives,  and  if  th»*  matrix 

r 

A  -  3U  A  Cu  -  EQ  i  -  (tr  *)1  *  B:  8*'(0)CI;  where  ?*'(.0)  denotes  the 

Jacobian  matrix  of  the  mapping  7*  evaluated  at  the  origin,  is  positive 
definite,  then  the  solution  x  =  0  of  Equation  (J*-6)  in  unstable. 

Proof-  The  proof  proceeds  ir.  the  same  manner  as  that  of  Corollary  J. 
Q.E.D. 

We  now  give  several  examples  of  the  application  of  the  above  results. 


}.  Example  1 

For  our  first  example  we  consider  the  distributed  network  of  Figure 
.  .2  which  was  examined  from  the  viewpoint  of  linear  network  theory  in 
Chapter  I.  The  network  is  redrawn  in  Figure  b . ]  to  show  explicitly 
that  it  is  a  member  of  the  class  of  networks  having  the  form  of  Figure  U.j. 
Uote  that  we  have  replaced  fin-  resistor  r  which  was  characterized  by 
the  equation  i  *  gv  by  a  lumped  memory less  nonlinear  clement  character  ...  ed 
sy  the  equation  1  -  f (v) .  V»  assume  that  the  fur.  ’tier;  f  satisfies  a 
Lipschitz  condition  in  some  neighborhood  cf  v  =  0,  and  that  i‘(0)  •  0. 

Tnas,  the  resistor  r  of  Chapter  1  is  but  a  special  case  of  the  type  of 
*■  ieinent  we  shall  consider.  -.  r  th.s  network  Equations  (E-l)  and  (l»-2) 

-  -or.o: 
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Figure  5.1.  Network  for  Example  1. 


[C ]x^( t)  ■  [0]x1(t)  ♦  [-1 


0] 


iBU) 


VQ  (t) 

1 

1 

O 

O 

O 

‘a,(t)  I 

i 

1 

V° 

m 

1 

x:(t)  ♦ 

0  0  0 

1,  (t) 

1 

1 

i  (t) 

0 

0  0  1/R 

v_  (t) 

\ 

U  m 

W  mi 

as,  P  -  [C],  A  -  [0],  B»  -  [-1],  B  „■  (l  0),  8.  ■  [l],  B  - 

t  i 

-  tl],  C  -  [1  or,  c6  -  Ui,  C  -  [OJ,  D:  »  to],  d6  -  to], 

_!  1 

=  [i'B].  Also,  B_?*(C.x)  *  -l'(x.),  fl-  and 


A'BuACn  +  EutEn*  ^ f  r  ^  > ‘ n  f-r 


5|H  - 


VrrfTTT,  .  yj\ 


Thus,  by  Theorem  1  of  Section  the  solution  x.  *  i  T.-y>»j.»nt. icel  l.y 


stable  if  there  exists  K  >  0  such  *h-«t 


j  / CR  ♦  20)  -  J | K  -  2  It 

V  r  r-> :  20i  /*:  *  s(*i))>0> whcn  i\i  - *■ 


In  a  similar  manner,  we  may  apply  Theorem  2  of  Section  •>-?.  to  obtain: 
The  origin  is  unstable  if  there  exlstn  T>  0  such  that 

j  |(R  *  Z0)  ♦  |R  -  2  | | 

‘i  f‘  r((R  ♦ . zo)  ~|r  - "zo|  Jxi "  f(xi:)  >  0>  when  i xii  *  r 

These  criteria  may  be  specif i*  1  graphically  as  in  Figure  If, 

in  some  neighborhood  of  the  origin,  the  function  !'  lie;;  within  the  open 


r>  ) 


S,  ii'nble  /  /' 


Vs/V/V/A 


'  /  / ,  / 
'  /  /  / 

'  /  s  / 


/////,  i> 


Vn> 


•  .*■  a  for  Exarr.;. 
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region  labeled  stable  in  igure  5.2  (as  curve  does),  then  'he 

equilibrium  solution  -  0  i  r.  asymptotically  stable,  if,  in  some 
neighborhood  of  the  origin,  f  lies  within  the  open  region  labeled  unstable 

(■mrye  0,  for  example),  the  solution  x J  *  0  la  unstable.  If  the 

function  f  lies  within  Uw  remaining  region,  *he  solution  x^  »  O  may  or 
may  not  be  stable.  Also,  according  to  Corollary  i  of  Section  4-2,  if 
for  all  x^,  f(x  )  lies  Within  the  open  region  labeled  stable,  then  the 

solution  x^  ■  0  is  completely  stable.  It  was  mentioned  in  Chapter  I  that 

results  of  this  type  would  be  obtained. 

There  are  two  cases  to  consider  in  determining  the  stability  and 


instability  regions  of  Figure  5.2.  If  R  >  Zq,  then  the  straight  lines 
L ^  and  l ^  are  determined  by  the  equations 


As  mentioned  In  Chapter  I,  we  note  that  as  2^-yR  the  lines  i ,  and  i,} 
both  approach  the  line  x ( x 1 )  *  (-l/R)x^.  This  results  in  stability 
and  Instability  reg'or.s  whi-h  approach  those  shown  in  Figure  l.^a-  One 
would  expect  to  obtain  s.,-h  regions  for  the  lumped  network  of  Figure  5-2, 
having  no  transmission  iir.e.  Of  course,  when  Z  *  H  the  network,  of 

V 

Fi.r-i-";  5-2  is  equivalent  to  s.i-h  a  lumped  network. 
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4.  Example  g 


In  this  example  we  consider  a  network  which  contains  two  memory leas 
nonlinear  elements.  The  network,  shown  in  Figure  5.3,  consists  or  two 
lumped  network*  connected  together  by  a  lossless  transmission  line.  The 
lumped  network*  are  identical  except  that  the  nonlinear  elements  are 
characterized  by  (in  general)  different  functions  f  and  f^,  and  the 
values  of  and  need  not  be  equal.  Me  eesume  that  and  satisfy 
Lipschitz  conditions  in  some  neighborhood  of  the  origin,  and  that  f^(0)  = 
f^(0)  ■  0.  For  this  network  Equations  (4-1)  and  (4-2)  become: 


Figure  5- 3’  Network  for  Kxample  2. 
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ftms,  by  Theorem  1  of  Section  5-2,  the  volution  x  «*  0 
is  asymptotically  stable  if  there  exists  sene  K  >  0  such  that: 
When  Yo  >  g, 

xj  “  g+T*  1  ~"2g"g 1*  *J  +  f^xj^  >  0  '  for  txijl<K>  3ml>2  * 


when  Yo  <  g  , 


xj  ■  |(*-)3  xj +  f^x.5^ >  °  *  for  ixji<*  >  • 


Similarly, 


A  -  Bl  A  Cd  -Enf  E“  -  (tr  f  )Ifl  =  A  +  G  A  -  3T 


2  ,  „2 

l-G  ♦  {  *  J* 


iv«LrV! 


V* 


-G  ♦ 


J_  i  c^|  ly«l  *  vm 
g+Yo  “  4  g  \  Yo*g  |  J 


-G  + 


g+Y, 


2g-3Y 

2(;_ 


< 


-G  - 


1  G! 


2  l  g+Y, 


0 

°  g+Y0  \  2g  |  J 


,  if  Yo  >  g  , 


i  g! 


'  J  '  2U+Yo/J 


,  if  Y0  <  g  . 


Thus,  by  Theorem  2  of  Section  ‘>-2,  the  solution  x  *  5  is  unstable  if 
the r  *  exists  same  y  >  t  su-'h  thxt: 
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♦ 


When  Yo  >  g, 

xjito  ‘  &&  ]  xj  *  W}  <  0  fcr  ^<r  '  *ml’2  ' 

When  *o  <  g , 

V10  * « xi  *  W <  0  ■ for  l\il<r'  •)*1’2  • 

As  in  the  first  example,  the  above  criteria  may  be  specified  graphic¬ 
ally.  We  may  draw  a  figure,  identical  to  Figure  5.2,  and  require  the 
curves  of  both  f.  and  fg  to  lie  within  the  open  region  labeled  stable 
to  ensure  complete  stability  of  the  solution  x  =  5.  Similarly,  if 
in  some  neighborhood  of  the  origin  the  curves  of  both  and  fg  lie 
within  the  region  labeled  unstable  then  the  solution  x  »  0  is  un¬ 
stable.  The  lines  and  ig  are  now  determined  by  the  following 
equations: 


If  Yq  >  g,  then 


il(xj) 


[G 


_G__ 


2g+5Y 

^  "  2g  ^]xj 


-  [G  - 


2g-5Y 

2g 


o 


If  Yq  <  g,  then 

W  *  •  [0  •  I  (ife;)]x3  - 

12(,3)  •  -  lo  *  2  ■ 

We  note  that  as  r  -«  0  the  lumped  networks  tend  to  become  uncoupled 
from  the  transmission  line,  and  hence  from  each  other.  Also, 


3 


T9 


P-*  0-e  g-»  •»,  and  each  of  the  above  Haas  approaches  the  line 
“V  ■  -  °Y  8uch  a  result  is  satisfying  since  when  r  ■  0 
our  network  becomes  two  staple  networks  of  a  type  already  considered 
for  which,  f(xj )  >  -  Gx^  -«♦  complete  stability  and,  f(x^)  <  -  Stx^  “> 
instability. 

If,  in  this  example,  the  two  nonlinear  elements  were  not 
independent  of  one  another,  it  might  be  more  convenient  to  use 
Corollaries  2  and  3  of  Section  5-2  to  determine  asymptotic 
stability  or  instability  of  the  solution  x  »  0.  Suppose  the 
nonlinear  multiport  of  Figure  5.^b  is  characterized  by  the  equa¬ 
tion 


where  5  (5)  *•  5,  and  $  satisfies  a  Lipschitz  condition  In  same 
neighborhood  of  the  origin.  Since  0  we  have  a  s  3  .  It 

would,  in  general,  be  rather  awkward  to  try  to  obtain  graphical 
stability  criteria  as  before.  If,  however,  for  all  x  in  seme 
open  set  containing  the  origin,  -  {x:x  £  En,  |fx||  <  K),  the 
function  8  (Cj  x)  has  continuous  second  partial  derivatives 

and,  if  Q  '(5)  denotes  the  Jacctian  matrix  of  the  mapping 
evaluated  at  the  origin.  then,  according  to  Corollaries  2  and  3 
oi  G*. ction  5-2:  If  the  matrix 


..  p....itive  definite,  U*«  solution  *  •  8  it  asymptotically  stable.  If 
».hp  matrix 


r  j£_  3  | 

*G  *3r  ■  t  •  t  r 


V«l  ♦  V* 


Y  ♦£ 


1 


5  '(<5) 


O  “  »  o 

is  positive  definite,  the  solution  it  •  6  is  unstable.  Such  criteria  as 
these,  even  for  a  much  larger  network,  should  not  be  too  difficult  to 
verify,  provided  that  adequate  computing  facilities  are  available. 


Example  3. 

In  this  example  we  consider  a  large  network  consisting  of  an 
arbitrary  (finite)  number  of  voltage  controlled  nonlinear  resistors, 
bavin.;  capacitance  in  parallel,  connected  together  in  an  arbitrary 
manner  by  lossless  transmission  lines  having  lumped  resistance  at  each 
end.  The  elements  of  the  network  are  shown  in  Figure  5.U.  We  assume 
tha*  all  transmission  lines  in  the  entire  distributed  network  have 
tr.-j?  i'ame  characteristic  impedance  Zq.  The  parameters  a  and  £,  however, 
may  pe  different  for  different  lines.  We  let  n  denote  the  number  of 


(a)  (!,. 


(a)  Typical  jumped  netwoiv.  (b)  Typi.-al 
■  ■  ••rcc.uv.-ct ing  line 
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Figure 


lumped  network*  in  our  large  <Jl<n.rlbuted  network,  and  usr  the  subecript 
j,  ... t  n,  to  denote  each  particular  luaped  network.  Ve  denote 

by  the  positive  integer  the  number  of  Interconnecting  lines  which 
are  connected  to  the  j-th  luaped  network.  £n  Figure  5.5  we  show  a 
typical  l\xped  network  and  the  lines  connected  to  it. 


Figure  5.5.  A  r /pi cal  subnetwork  for  Example  3. 

The  Kirchoff's  voltage  and  current  law  equations  for  this  subnetwork 
are: 
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For  this  example,  the  limped  linear  multiport  of  Figure  k,2  con¬ 


sists  of  the  collection  of  all  lumped  linear  subnetworks  of  the 
type  shewn  in  Figure  ;>■!?•  Thus,  for  the  entire  network  vc  have: 


Equation*  (5-7)  are  not  necessarily  in  the  $a:  •  form  as  Equations 
(--1)  and  (4-2).  This  la  be--ause  the  elements  of  t;»  sectors  u  and  w 
may  not  be  in  the  proper  ord^-r  as  defined  in  Chepte r  tV  Let  us  investi¬ 
gate  this  matter:  Since  each  transmission  line  ir.  tlv  distributed  network 
is  connected  to  a  port  of  ih  lumpeu  linear  multip-  V  •'  r  which  the  port 
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volume  ha*  been  chosen  ar,  » hf-  independent  variable,  the  vector*  u 


and  v  of  Equations  (h-I)  and  must,  according  to  Chapter  IV, 

be  of  the  form 


U  a 

V 

i  * 

»  W  - 

kJ 

Im 

It  appear*,  therefore,  that  we  need  only  relabel  the  port  variable*  to 
-#  -* 

give  u  and  w  the  proper  form.  Thing*  are  not  thia  simple,  however, 

since  thi*  relabeling  cannot  always  be  a.: coup 11 shed .  We  see  that 

_« 

the  voltage  variable*  in  the  vector  u  occur  in  groups  such  that  all 
voltage  variables  associated  with  any  lumped  subnetwork  (e.g.,  the 
J-th  subnetwork  of  Figure  5.5)  appear  in  adjacent  locations  in  the 
vector  u  .  Hence,  if  these  variables  could  be  relabeled  as  required 
then  at  moat  one  subnetwork  would  have  so ne  ol  its  independent  port 
voltage  variables  contained  in  both  of  the  vectors  and  v^.  Con¬ 
sider  now  the  following  counterexample:  suppose  there  are  four 
lumped  subnetworks  and  six  interoonnec *ing  lines.  Let  the  inter¬ 
connections  be  made  such  that  each  subnetwork  ia  connected  to  one  end 
of  three  different  lines,  each  of  which  is  connected  to  each  of  the 
other  three  subnetworks,  as  indicated  in  Figure  5.6.  If  at  most  one 


:.V» work  :'.r  1  junterexamj. ;e . 
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Figure  ‘,.6. 


subnetwork  ha*  some  of  Its  independent  port  voltages  in  both  of  the 
vectors  end  v^,  then  the  port  voltages  associated  with  at  least 
two  other  subnetworks  are  all  members  of  either  or  vft.  However, 
since  the  port  voltages  associated  with  each  end  of  any  line  must  be 
labeled,  one  with  an  €  subscript,  and  the  other  witn  a  S  subscript, 
it  is  clear  that  if  any  two  subnetworks  are  considered,  all  six  ox' 
their  port  voltage  variables  cannot  be  relabeled  with  the  c  sub¬ 
script.  Similarly,  they  cannot  all  be  relabeled  with  the  0  subscript. 
Arriving  at  this  contridiction  proves  that  the  proper  relabeling 
cannot  be  accomplished  ior  this  network. 

In  order  to  put  Equations  (5-7)  in  the  form  of  Equations 
(4-1)  and  (4-2)  one  must,  in  general,  rearrange  the  port  voltage 
variables  in  the  vector  u*,  and  similarly  rearrange  the  elements  of 
w  .  If,  for  each  line  in  the  distributed  network,  one  end  is  chosen 
arbitrarily  and  the  corresponding  port  variables  relabeled  with  an 
e  subscript,  and  if  at  all  other  ports  to  which  a  transmission  line 
is  connected  the  port  variables  are  r*  ibeled  with  a  0  subscript, 
then  the  vector  u  may  be  put  in  the  required  form  by  simply  rearrang¬ 
ing  several  pairs  of  its  elements.  Performing  the  same  operation  on 
the  vector  w  will  put  it.  in  the  required  form,  let  U  denote  a 
matrix  which  performs  this  rearrangement  operation  on  the  vectors 
u  and  w  .  Ther. 

u  Uu  ,  and  w  ■  Uv  .  (5-8) 

V'c:  shall  now  briefly  inve  stigate  trie  nature  of  the  matrix  U. 
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If  one  starts  with  the  identity  matrix  and  interchanges  the 
i-th  and  the  J-th  columns,  the  resulting  matrix  Ujj,  when  pre Multi¬ 
plying  any  vector,  will  give  the  original  vector  with  the  i-th  and 
j-th  elements  interchanged.  In  the  same  manner,  we  may  construct 
a  nratrix  U  which  performs  the  operation  of  interchanging  several 
pairs  of  elements,  by  interchanging  the  corresponding  columns  of  the 
identity  matrix.  It  is  thus  clear  that  U  is  nonaingul&r  and,  more¬ 
over,  is  its  own  inverse,  since  UUx  -  Xx  for  all  x,  He  also  note 
that  if  A  is  a  diagonal  matrix,  then  UAU  *  A. 

If  in  Equation  (5-7)  we  define 
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I 


and  if  in  Equation  {5-6) 


I  0 


0  u. 


„ 

then,  since  u  *  Uu  and  w  ®  Uw,  Equations  (5-7)  become 

Pl(t)  ■  Ax(t)  ♦  B#Ud(t) 


Uw( t)  =  C*x(t)  *  B#Uu(t)  , 


#  * 

and  hence  (noting  tJD  U  »  D  )  , 


pj(t)  =  Ax(t)  ♦  B#Uu(t) 


v(t)  *  UC*x(t)  +  D*u(t)  . 


These  equations  are  in  the  form  of  Equations  (4-1)  and  (4-2)  if  we 


define 


-  r  *  1  *  , 

B  =  B  U  =  t  Bj  |  B  nUjj  ] 


C  =  UC*  =  - 

Va 


Now,  letting  G  *  1/H  and  e  =  jj  +  p  , 


=  g+Yo  » 


1**0? 


(g*Y0)‘ 
(y  *<o‘ 

u 


i 


for  all  l  (.  U 


AIho,  “  l|C^||  -  0  for  all  {£  Q  and,  hence 
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- 1 

» 

X 

g 

so 

1 

O 
_ _ _ 1 

•*  -x,o2 

\o 

O  V2_ 

p\ 

where  X  «  (g+Yo'"x 
Similarly#  since 


]V8[_ 


•  G  ♦  Y. 


0  (8ty 


•  G 


i  -  (y«)2/<y g)2 


lyel  +  V« 

Yo+« 


for  all  tea. 


BnTE 


<  iygi  *  ye 

1  ye 


and 


<tr  Y)Iq 


iyei  +  *0  +  g 

Yo+* 


where  denotes  t!.<:  total  riur.be r  of  lines  in  the  ■entire 
r.etowrk.  Therefor-: . 


distributed 
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4  -  V  cti  *  En  *  h  +  (tr  f)In  s  A  f  * 


?  /  tVd  *  Y  ♦« 


Y  +• 


kj+SN 


'  o 


O  ■■ 


k  +2N. 
n  f 


klg2  gc1+afi)og  i  |YQ-gf  ♦  ye 


■xw  r  e*Y0  T 


-knG  +  ^  + 


yQ+g 


*>2  (VfVo2  I  ly«l  ♦  Vs 


Vs 


For  j  «  1,  . . n,  let 

,2 


*,0  (k  ♦  w.xr  I  It  -el  ♦  Y  +g 

■j#*!*r*“1T5 — 1  t.7*- 


Then,  if  Yq  >  c  , 


PJ  *  -KJ 


G 

g+Y^. 


l  + 


i)]-. 


g2y 


g(g+Y0) 


If  Y  < 


"i 


+  N 


/  S+Y 


(5- 


(5-: 


As  r  -*  0  note  that  f-*  »,  and  in  Equations  (5-9)  and  (5-10), 

T-.  _  -k<G  =  -Y.  ,/R.  Causing  r  to  become  small  has,  of  course,  the 
•fleet  of  tending  to  uncouple  the  lumped  circuits  from  one  another 


*  -  «c  -  V  *S  •  (tr,)In 

kjO2  (k  >02  |  |Y  -g|  ♦  If  « 

ki°  *  jir  - -V~  ra — 


O 


o 


‘kn°  *  g+Y 


*/  <kn^v° 


!Vg'  ♦  V« 

V* 


For  J  **  1,  . . n  ,  let 


k.Q2 

(k,  .  2N,)02  l|Y  -8|  ♦  Y 

’j  ‘  -kj°  *  sir  • 

'tg 

\  Yotg  / 

Then,  if  Y  >  g, 

u 

r  2 

Ao 

qj  ■  -‘jI°  -  3ra 

I1  -  d 

“  g(g+Y^  * 

(5-U) 


If  Yq  <  g  , 


»[•-*(  AH 


G 

g+Y, 


(5-12) 


In  Equations  (5-11)  and  (5-12)  it  la 

\l  "* 

Sim.. 

BI  =  BI  U*^CI  * 

l...  application  of  Theorem;:  i  and  2  of  Section  5-2  yield  the  following 

stability  criteria:  The  solution  x  »  5  is  asymptotically  stable  if 

there  exists  sane  K>  0  such  that  x.(-p.x  ♦  f  (x  ))  >  6  ,  for 

J  J  J  J  .J 


also  true  that  as  r  -*  0, 


flixl^ 


f  (x  ) 

n  n ' 
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Jx.j  <  K,  j  =  1,  . ..,  n,  where  p.  it  given  fc>  Equation  (5-9)  for 

J  v 

Y  >  g,  and  by  Equation  (‘>-10)  for  Y  <  g.  'IYic  solution  x  ■  5  ii 
o  —  o 

unstable  if  there  exists  sane  )  >  0  such  that.  *j{-q,Xj  *■  ^ 

for  |*j {  <  r -  3  ~  1 ■»  -  “•  R»  “here  q^  is  given  by  equation  (?-ll) 
for  Y  >  g,  and  by  Equation  (5-12)  for  Y  <  g.  A*  in  the  other 
examples  we  may  specify  these  criteria  graphically.  For  J  ■  1,  ... 
we  may  draw,  as  in  Figure  5.7, 


J 

p.x 


J  J 


Figure  5-7*  Stability  criteria  for  Example  3. 


regions  in  the  x.,f.(x.)  plane,  ir  f  lies  in  the  open  region 
J  J  J  J 

labeled  stable  for  all  x,  /  0  in  some  neighborhood  of  the  origin, 

J 

the  solution  x  -  0  is  asymptotically  stable.  If  for  all  jrf  0, 
f  (x  )  lies  in  the  open  region  labeled  stable,  then  the  solution 

J  J 

x  a  0  is  completely  stable.  Similarly,  if  fj  lies  in  the  open 
region  labeled  unstable  for  all  x.  /  0  in  some  neighborhood  of 
the  origin,  the  solution  x  -  0  is  unstable. 


We  hav-1  noted  that  as  r  —0,  each  of  the  p^,  approach 
-*/R.  This  ir  oxa.-tiy  what  we  would  hope  to  obtain  since  when 


0, 


,  n* 
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r  »  0,  the  n  lumped  networks  are  uncoupled  end  they  ere  completely 
•table  If  fj(x^)  >  “^jG-  our  th*°**y  **  have  established 

criteria  which  ensure  complete  stability  (end  instability)  even 
when  the  lumped  networks  ere  coupled.  Rote  that  the  criteria 
art  Independent  of  a  end  I  for  every  line,  and  arc  also  Indepen¬ 
dent.  of  the  value  of  for  J  -  1,  ...»  n. 


Chapter  VI 
CONCLUSIONS 


A  stability  theory  for  nonlinear  distributed  networks  has  been 
presented,  in  the  development  of  this  theory  three  major  steps  were 
taken,  ttrst,  it  was  shown  that  the  electrical  behavior  at  the  ports 
of  certain  two-port  networks  containing  lossless  transmission  lines 
may  be  described  by  a  system  of  functional  equations.  Next,  a  class 
of  nonlinear  distributed  networks  vna  defined  and  it  vhs  shown,  using 
state  variable  techniques,  that  the  behavior  of  any  network  in  this 
class  n»y  be  characterized  by  a  system  of  functional -differential 
equations.  Finally,  a  Lyapunov  functional  was  presented  and  the 
stability  theory  for  functional -differential  equations  was  used  to 
obtain  several  theorems  and  corollaries  which  specify  sufficient 
conditions  to  ensure  that  the  equilibrium  state  of  a  given  network 
in  the  above  class  is  stable,  asymptotically  stable,  completely 
stable,  or  unstable. 

It  has  been  shown  that  the  stability  criteria  that  this  theory 
may  specify  for  any  particular  distributed  network  are  independent  of 
the  length  of  the  transmission  lines  and  also  independent  of  the  phase 
velocity  (  a  l/a)  of  the  lines.  Generally,  the  criteria  are  also 
independent  of  the  values  of  the  reactive  elements  contained  in  the 
network. 

The  stability  criteria  nay  not  always  be  the  best  that  one 
might  hope  to  obtain.  For  example,  in  the  paper  by  Brayton  and 
Mlranker  [10],  which  is,  at  this  time,  the  only  other  comparable 
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tt-.'  ify  known  to  the  author,  -i  i  ibUlty  criterion  for  Uv'  network 
of  Example  1,  Chapter  V,  is  given. 

They  obtain:  If 


t  >  - 


$  it 

K  * R 


(6*1) 


then  the  aolution  ^  o  is  completely. atahle.  Ijy  application  of 
Corollary  1  to  this  example,  wo  obtains  It  z  -£  R  and 


f  > 


P.  •  2C  ’ 


or  if  2  5  R  and 


f*  >  "  -  ~ 
72  K 
Zo  ° 


(6-2) 


(6-5) 


then  the  solution  0  la  completely  stable.  It  Is  easily  shown 
that  if  the  ratio  f/2q  lc  leon  than  approximately  0 .648  then 
l./R  -  2/Zo  >  -  R/(Z2  +  Ii2)<  and  hence  the  "tit.erion  of  Equation  (6-1) 
is.  less  reatrictlve  (in  terns  or  admissible  functions  f)  than  the 
criterion  of  Equation  (6-2).  in  cane  the  ratio  p/2  is  greater 

r,'<  ■  •*!  'y, 

than  approximately  j.  34. then  R/Z*"  ■  2/2  >  -  R/(  7?  *■  F‘ '  .•  and  hence 

O  O  C/ 

the  criterion  of  Equation  (6-1)  is  leas  restrictive  than  the 
criterion  of  Equation  (6-5)..  if,  however 4 
0.648  <  r/zo  <  •!  .'A, 

then  the  criteria  of  Equation:;  it—;.*)  end  (6-5)  ate  less:  restrictive 
t nun  the  criterion  of  Equation  1 '•-!).  For  example ,  if  7.  10  it 

an  i  H  *  9  At  then  our  criteria  imply  complete  stability  ;f  t‘  > 

•;.()  :Vj-,  Brayton  and  Mi '-.'inker'  .  criterion  Implies  omj>!  -.■*•■  stability 
)'  i  *  >  -  0 .0497  • 
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Several  r  a.;  f  future  research  along  toe  liner;  of  this  work 
prove  fruitful.  Kir  .a  of  all,  it  wculu  be  significant  to  extend 
t-'.e  of  iittrieutva  networks  which  say  he  described  by  functional* 

differential  equations  (and  hence,  to  which  our  stability  theory  night 
apply)  beyon-i  the  class  which  we  have  defined.  One  such  extension 
•sight  proceed  alone  the  lines  of  removing  the  restriction  in  condition 
concerning  the  zeros  of  transmission  at  s  «  ■*  of  the  matrix  D. 

If  thin  could  fce  ione,  th>-n  networks  containing  transmission  linos 
terminated  at  one  end  by  only  lumped  memoryless  nonlinear  elements 
could  be  studied.  We  'Quid  study,  for  example,  the  network  of 
figure  6.1.  The  stability  of  the  equilibrium  state  v  ■  0  for  this 
network  cannot  be  studied  by  our  theory. 


Figure  o.l.  Another  nonlinear  distributed  network. 


Another  extension  of  our  work  might  be  to  consider 

networks  contain: my  types  of  distributed  elements  which  are  more 
r.'.'ral  than  1<. ■  L";  ’■  runsmisa ion  lines .  A  theory  for  networks 
■on tain l ng  P.C  i i n-;  • ,  l'  -r  example ,  should  be  quite  uieful. 

The  hy-i}:;..’.  t ....  cf  Equation  {'  *1)  i-  probably  not  the 

..'-y  >ne  ;h  ~  •  !  :  f  .no  that  lea  is  to  a  useful  stability  theory. 
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Of  course,  the  conception  of  Lyapunov  functionals  is  of  an  art 

than  a  science.  If,  however,  one  is  clever  enough  to  fina  ther 
Lyapunov  functionals,  it  shoui d  he  rather  routine  to  develop  new, 
and  perhaps  better,  stability  .-.riteria  for  the  class  of  distributed 
netverks  which  ve  have  defined. 

Networks  containing  only  lumped  meoorylese  elements  and  t  ma¬ 
nias  ion  lines  (i.e.,  no  reactances)  can  certainly  also  have  stable 
and  unstable  equilibria.  A  stability  theory  for  such  networks  should 
be  a  valuable  contribution. 

Finally,  having  defined  certain  two -ports  by  functional  equations, 
one  is  lead  to  conjecture  about  the  perhaps  academic  problem  of  devel¬ 
oping  a  theory  of  analysis  of  networks  containing  what  might  be  called 
"functional  elements".  A  functional  one-port  might  be  described  by 
an  equation  of  the  type 

i(t)  --  f(vt), 

where  i(t)  denotes  the  value  of  the  current  through  the  element  at 
the  time  t,  and  v  is  a  function  (the  voltage  across  the  element) 
on  some  time  interval  for  which  t  is  the  right-hand  end.  Perhaps 
the  theory  of  dynamical  systems  would  be  the  proper  setting  for  such 
a  problem.  Choosing  the  proper  state  space,  and  defining  the  order 
o*’  .omplexity  for  such  networks  floes  not,  at  this  time,  seem  to  be 
a  trivial  problem. 
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APPENDIX  A 


In  this  appendix  we  study  some  of  the  aspects  of  the  problem  of 
locating  the  zeros  of  the  function 


Y(s)  a  +  sC  +  =- 
r  Zo 


(ft  +  Zc)eTS  ‘  (H  -  Z0)e~** 

(R  ♦  Z.  )eTB  +  (R  -  Z0)e‘TS 


The  equation  Y(s)  >0  is  the  characteristic  equation  of  the  distributed 
network  considered  as  an  example  in  Chapter  I.  It  is  clear  that 
Y(s)  a  0  if  and  only  if 

(ZQ  +  sCrZQ)t(ft  +  Z0)e2ts  +  <R  -  ZQ ))  «  -r(R  +  Zq)c2ts  +  r(R  -  ZQ)> 
or 

t(r  +  Zq)(R  +  ZQ)  +  sCrZo(R  ♦  Z£J)]e2ts  *  t(R  -  ZQ)(r  -  ZQ)  -  sCrZ0(R  -  ZQ)J. 

Let  ua  define  z  ■  2%&,  and  let 

CrZ0(R  ♦  Z  ) 

Q  =  - - -  , 

3  *  (r  +  ZQ)(R  +  Z0)  , 

CrZ0(R  -  Z0) 

I  a  *  ■  i 

2t 

6  *  (P  -  20)(r  -  ZQ). 

It  follows  that  Y( s )  =  0  if  and  only  if  s  =  £/£ . ,  where  £  is  a  root  of 
the  equation 

(az  +  0)eZ  -  yz  +  6  •  (A-l) 

The  location  c-f  the  roots  of  Equation  (A-I  )  has  been  studied  exten¬ 
sively  by  E.M.  Wright  We  shall  adopt  some  of  his  techniques 
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here.  We  first  trams form  Equation  (A-i )  into  a  simpler  equation  by 
replacing  x  by  a  -  la(a/rj.  Rote  that  a/r  ■  -(R  +  ZC)/(R  -  zc)  ■ 
-1/r,  where  T  is  the  reflection  coefficient  at  the  right-hand  end  of 
the  transmission  line  in  Figure  1.2.  In  case  a/r  >0  (l-e.,  r  <  0, 
which  is  the  case  for  the  example  of  Chapter  I)  we  obtain 


where 


(x  -  A  ♦  B)e*  -  x  -  A  -  B, 


(A -2) 


A  -  -  e/£L±-fiZr  +  laja/r|  -  -  g  +  ln|-l/r|, 

B  .  B  is-  . 

2  CZQ 

Suppose  r  =  -1.78*  R  =  1.42,  C  *  1,  Zo  «  2.19,  T  =  s/6.  It  then  follows 
that  ln|a/r|  *  1.55,  A  as  2.14,  and  B  »  0.48.  If  we  define  the  function 
Cl(B)  by  cx(B)  «  ln(B  +  1  +  V  B2  +  2b'  )  +  V  B2  +  2b\  then  c^O.48) 

2.05*  According  to  Wright  [5],  if  B  >  0  and  A  >  c^(B)  then  Equation 
(A -2)  has  exactly  two  real  roots.  Furthermore,  all  of  the  roots  of 
Equation  (A-2)  have  real  parts  which  are  less  than  or  equal  to  the 
value  of  the  larger  of  these  two  real  roots.  It  is  clear  that  this 
statement  also  applies  to  the  zeros  of  Y(s)  for  the  shove  parameter 
values . 

We  shall  now  show  that  the  two  real  zeros  of  Y(s)  lie  in  the 
left  half  of  the  s -plane  for  the  parameter  values  given  above.  This 
will  then  prove  that  Y(s)  has  no  zeros  with  positive  real  parts.  We 
first  rewrite  Y(s)  aa 


Y(s)  =  7  *  sc  +  !  Srr - -  • 

r  ZQ  i  e"To  +  r 
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1  i  f  e**  „  r  1 

Then,  if  we  define  y^s)  =  -  -  -  sC  and  y^s)  -  -r 7-  I  “jjjj - -  I  » 

it  follows  that  Y(s)  o  if  and  only  if  y,(s)  «  yD(o).  On  the  rani 


axis  in  the  s -plane  we  have  ^  yg(a)  *  [(a^  +  F)2t*‘ 

o 


2T8 


(e2*8  -  r)2re^a]  /  (ecTS  +  r)*  - 


2f  a 


c'TS 


4rter 


Z0(e2rs  +  r)2 


Thus,  r  <  0=*  y,  (s)  <0  for  all  real  s.  Since  we  also  have  y~(-«)  •« 

uS  d  • 

-  —  >  y2(o)  -  r(rrr) !  r'  y2(*°  55  \r> the  function 

"o  c  o  o 

behaves  as  shown  in  Figure  A-l,  for  real  s.  The  point  s  -  a  is  the 
solution  of  the  equation  e210  +  T  -  0.  The  important  thing  to 


Figure  A.l.  The  functions  y^s)  and  y2(s)  for  real  s,  when  T  <  0. 

note  here  is  that  y  ( s )  is  a  strictly  monotonic  decreasing  function 

<u 

sad  —  y  (s)  is  strictly  monotonic  increasing  for  s>  0}  hence,  the 
ds  J  2 

straight  line  y,  -  1/r  -  intersects  this  curve  in,  at  most, 

•,~o  points.  Tl.i^  agree.-:  with  Wright’s  results.  In  our  case  there 
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are  two  interne  t  ions  and  they  occur  at  approximately  a  »  -O.819  and 
•  «  -0.2k2.  Ttoat  is  ,  Y(s)  has  exactly  0*0  real  zeros  and  they  are 
both  negative.  Hence,  ve  have  shown  that  for  the  parameter  values: 


r  *  -1-76,  R  «  1.U2,  C  -  1,  ZQ  •  2.19,  and  t  =  x/6,  the  zeros  of  Y(a) 
all  have  negative  real  parts. 

The  aain  reason  that  Figure  A.l  has  been  introduced  is  that  it 
can  provide  a  certain  amount  of  insight  into  the  stability  problea  for 
the  linear  distributed  network  of  Figure  1.2.  We  shall  not  consider 
this  topic  in  great  depth,  however  the  following  remarks  seem  to  be 
appropriate.  First,  we  saw  in  Chapter  I  that  an  unstable  network 
was  made  stable  by  simply  increasing  the  value  of  the  capacitor  C- 
In  Figure  A.l  we  see  that  varying  the  value  of  C  simply  causes  the 
slope  of  the  straight  line  representing  y1(s)  to  vary.  Hence,  when 
l/ZQ  <  -l/r  <  l/R ,  we  can  always  make  C  small  enough  that  Y(e)  will 
have  two  positive  real  zeros.  Conversely,  as  we  increase  C  we  cause 
the  real  zeros  to  vanish  and  then  reappear  in  the  left  half  plane. 
Similarly,  we  can  cause  positive  real  zeros  to  exist  for  any  value 
of  C  >  0,  l/Z0  <  -  l/r  <  l/R,  by  simply  causing  the  point  a  to  move 
close  enough  to  the  e  =  0  axi3 .  This  corresponds  to  an  Increase  in 
the  value  of  t  (r  =  al),  which  in  turn  is  caused  easily  enough  by 
increasing  either  the  value  of  a  (a  •VvT)  or  f.  We  can  also 
force  the  real  ceros  to  occur  in  the  left  half  of  the  s -plane  by 
using  a  short,  enough  transmission  line  ( small  l)  cr  by  making  a 
small  enough. 

Let  us  now  consider  the  location  of  the  real  zeros  of  Y(s)  when 

-  l/r  >  1/R  (the  lumped  network  of  Figure  1.1  is  unstable  if  and  only 
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if  r  and  ft  have  values  such  that  this  inequality  is  satisfied,  since 
-  l/r  >  i/fc^s^g  <  -G).  If  T  <  O,  it  is  obvious,  from  Figure  A*l, 
that  Y(s)  will  always  have  one  positive  real  zero  when  -  l/r  >  l/R. 
If  T  >  0,  then  —  y  (a)  >  0  for  all  real  a,  and  hence  y(a)  behaves 
as  shown  in  Figure  A.  2.  Again,  we  see  that  Y(s)  will  always  have 
a  positive  real  zero  when  -  l/r  >  l/R.  Bence,  if  the  lumped  network 
of  Figure  1.1  is  unstable,  then  so  is  the  distributed  network  of 
Figure  1.2. 


Figure  A. 2.  The  functions  y1(s)  and  y2(s)  for  real  s,  when  F  >  0. 
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APPENDIX  b 


In  this  appendix  we  shall  prove  that  the  linear  space  1  m 
C(  ( -*,Q] ,  1°  )  with  the  compact  open  topology  is  metrlzable  with 
metric  p  defined  as  follows:  let  (t^)  be  •  sequence  of  real  numbers 
with  0  «  tQ  <  t^  <  .  . .  <  t^  <  t^+i  <  . . .  and  such  that  lim  t^  *  **. 

k-t  aj 

For  a  fixed  real  number  b,  0  <  b  <  1,  and  for  every  qj,  f  (.  C,  let 

«  k  » 

p(q>,$)  =  £  «*»  where  m^  -  min{  b  ,  sup(  |)q>(t)  -  i(t)||  s 

k-o 

-t^i  $  t  <  -t^  ]  ).  Since  0  $  ^  J  b^  for  all  k,  and  since  0  <  b  <  1, 

m 

it  is  clear  that  £  m.  always  converges  and  hence  p($>$)  is  well 
k-o 

defined.  We  first  verify  that  p  is,  indeed,  a  metric  on  C;  that  is 
+  9 t  U  I  £  C,  p  satisfies  the  three  properties: 

1)  o(?if)  2:  0,  p(p,$)  -  0  if  and  only  if  p  * 

2)  p($,t)  -  oihv)i 

5)  p($>f)  <  p(?A)  ♦  p(x,*). 

Properties  1)  and  2)  are,  of  course,  obvious.  To  prove  property  5) 
we  need  the  leana: 


If  A,  B,  C,  D  are  nonnegative  real  numbers  and  if  A  £  B  +  C, 
then  min(  A,D  ]  <  min{  B,D  }  +  roin{  C,D  }. 


proof:  If  D  <  B  and  D  <  C, 


raln(  A,D  )  <  D  <  D  +  D  =  mini 


If  D  <  B  and 
min{  A,D  )  g  D  $ 
If  D  >  B  and 

rr, i n {  A,D  )  <  D  S 


D  >  C, 

D  +  'T  -  m  I  :i  { 

D  <  ■, 

I;  ♦  t)  *  mir.  i 


B,D  ]  +  min(  C.D  )• 
« 

P,D  }  +  mint  C,D  }. 

B,D  )  ♦  mint  i:  ,D  j 


10b 


If  I)  >  I!  utld  U  >  ■"  i 


ttln(  A i0  )  <  A  <  B  ♦  C  -  ainf  B,D  )  ♦  mini  C,D  }  . 

Having  exhausted  all  pon s  1 t> i  li  ties ,  the  learn  is  proved.  Q.f.0. 

Nov,  for  any  k,  if  c  is  an  arbitrary  point  in  the  interval 
then 

Il9(tf)  -  <(o)||  -  t(9(o)  *  i(o)  +  X(o)  -  *(dl||  <  ((9(0)  -  x(c)l 

+  (|x(o)  -  *U)II  <  ^up(  j|cp( t )  -  X(t)||:  -tk+1  <  t  $  -  tk  )  ♦ 

oup(  |(x(t)  -  ♦(t.)||  :  -tk+i  <  t  $  -tk  )  . 

But,  since  a  is  an  arbitrary  point  in  [-tk+,,  -tjJ,  ve  therefor*  have 

sup{  |)jp(t)  -  i(t)|(  :  -tR+1  <  t  <  -tk  )  <  aup(  (J^(t)  -  x(t)||: 

*tk+1  <  t  <  -tR  }  +  sup(  ||x(t)  -  *(t)||:  -tk+1  $  t  <  -tk  ). 

By  the  above  lemma ,  ve  then  have 

min{  bk,  sup{  j|^( t )  -  J(t)||  s  -t^+.  S  t  <  -tk  )  )  $  aln(  bk,  aup(  |9<t) 

-  X(t)||  :  -tk+1  5  t  <  -tkJ  }  +  ain[  bk,  sup{  (|x(t)  -  ♦(t)j|:  -t.^  5  t 

<  -tk  }  },  for  all  k.  Thus,  <  p<9,X)  +  p(Xif).  Q.E.D. 

Let  3c  denote  the  compact  open  topology  for  C.  The  definition 

of  3c  is  given  in  the  first  footnote  of  Chapter  II.  Then  (C,  8C) 

denotes  the  topological  space  consisting  of  C  vith  the  compact  open 

topology.  If  p  is  the  metric  on  C,  defined  above  (  for  any  fixed 

b,  0  <  b  <  1),  let  m  denote  the  metric  topology  on  C.  Then  (C,  «,  ) 

0  p 

ui -notes  the  topological  space  consisting  of  C  vith  the  metric  topology 
"t  If  9  C  Ci  I'-'t  si.p.u)  l  denote  the  open  u-sphere  about  9, 

•  e • ,  S(9,u)  *  i  o  ~'tp,Y)  <  v  }. 

.'..corem.  The  ii.jM.iug: .- 1  .jiuce  (C,  3„)  is  met ri sable  vith  the  metric 
,  defined  above 
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proof:  Let  denote  the  base  for  9C  cone la ting  of  finite  inter¬ 
sections  of  members  of  (1  (see  first  footnote  of  Chapter  II).  Let-^ 
denote  the  bead  for  ^  consisting  of  *11  open  ^-spheres  in  (C,  i\p). 

We  prove  that  (C,9e)  *  (C,si^)  by  showing  that  if  $  is  a  point  in 
an  element  Bfi  of-^,,  then  there  is  on  element  Bp  of  containing 

the  point  9  and  contained  in  E,  and  conversely. 

a 

Let  B  be  a  member  of  4 I  .  Thus,  B„  *  ft  A(K.  ,  U.  ),  where 
c  c  c  W1  1  1 

A(Ki,U1)  la  the  aet  of  all  9  £  C  which  imp  K^,  a  compact  subset  of 
(-•,0],  Into  >  an  open  subset  of  Ea.  Let  9  €  Bc.  We  will  show 
that  there  exists  an  open  p -sphere  about  9,  S(<p,p),  which  is  contained 
in  Bc.  For  i  «  1,  . m,  is  compact  in  therefore  3  an 

integer  k  >  1  such  that  -tk#  <  t  for  every  t  €  K^.  Also,  9[KJ  is 

n  »  v* 

compact  in  E  .  For  every  y  £  9[K.  }  3  6  ,  0  <  t  <  b  ,  such  that 

1  Jr  j 

||y  *  -  y ||  <  .  The  family  of  6y/2  neighborhoods 

{  Ng  /0(y)  :  y  £  )  1b  an  open  cover  of  9(^1  and,  since  $[K«] 

y/ 

1  j  compact,  3  a  finite  subcover  (  N&  ^(y^),  •••  •  f)  ) 

yl  v<  " 

of  Note  that  ($j)  c  Ut  for  J  •  1,  ...  ,  f.  Let  &i  - 

yJ 

«in(  6..  /2,  ...  ,  5V  /£  ).  Then,  If  t  £  S(«p,  5.)  we  have  that,  for 

*1  *  g  «  1 

K  t  1  0,  ...,  k*-l  )  :  min(  bk  ,  sup(  ||ip(t)  -  *(<■)!' : )  } 


k  -1 


£  min{  bR  ,  sup(  ||qp( t )  -  v(t)j,:  -tR+1  <  t  <  )  )  < 


X-0 


k*-l 


nr.ir. (  bk,  sup(  ||<p( t )  -  5(t)|!.  S  t  <  -tR  )  )  ^  tr.k  < 


\<yt> 


r-’  „  u* 

^  \  ■  p(<p>7)  <  •  Th.-r.  fore,  v  t  <  K(  ,  min(  b  ,|(f(t)  -  f(t)||  )  < 

k*o 

6j  •  &ut  St  <  bk  {^(t)  -  i(t)||  <  6t  y  t  t  Thus,  ||f(t)  -  y^|j 

=  }}l{t)  -  p(t)  ♦  <p{t)  -  y.j;  5  -  p(t)||  t  |[i(t)  -  y^H  <  ♦  6y  /2 

£  5  /£  +  ft  /?  -  5  “>  *(*  )  €  K  (y.)  c;  U  .  ?tov,  let  p  - 

yJ  yJ  yJ  \  J 

min{  6, ,  ...  ,  t>m  )}  then,  S(qj,p)  1b  the  required  p-sphere,  since 

If  i  £  S(p,p)  and  If  t  t  Kj  for  any  it  (i,  ...  ,  m)  then  ♦(t)  £ 

that  is  f  maps  into  ,  thus  }  t  Bfi. 

Let  Bp  be  a  member  of©.  Thus,  *  S(p*,p*)  for  some  e  C 

and  some  p*  >  0.  Let  q>  c  B^.  Since  8(<p,p*)  is  open  in  (0,1^),  3 

p  >  0  such  that  S(^,p)  c  3(p',p')i  therefore,  it  is  sufficient  to 

m 

show  that  3  a  set  of  the  form  fl  AOC^,  Ui)  which  contains  9  and  is 

i=i  m 

contained  in  S(<p,p)-  In  case  p  >  Z  b  ,  we  may  simply  let 

k*o 

H  A(Kt,  Ut)  -  A(KX,  l^),  Where  1^  =  (t:  -1  <  t  $  0  },  -  tP. 

A(Ki»Ui)  =  C  and  hence  contains  p.  Also,  if  f  c  ACK^,  li^)  then 

p(ip,f)  <  2  bk  <  p,  hence  if  1.  S(^  ,p ) ,  «=<*  A(K,  ,  U,  }  c  S(p,p).  If 

k-o  A 

®  k  *  k* 

p  £  Zb,  then  let.  the  integer  k  >  1  be  chosen  ouch  that  b  <£ 

k=o 

u(  1  -i>)/2.  Consider  the  close  a  interval  [ -t^#,  Oj.  Since  p  is 

continuous  on  [  -t^*  ,  Oj,  V  t  i_  0]  3  6^  >  0  such  that  if 

-t,  *  <  t  ’  <  0  and  1 1  ’  -  1 1  <  6t  then  ||p(  t ' )  -  p(  t )  j|  <  .  Let 

"  i*k 

K.  (t)  denote  t.l.e  ..per.  6,  neighborhood  of  t.  Clearly,  N-  c 

t  °t 

(  t  '  :  |t t  |  <  b.  )  ■  Tf  f  i.~.i  iy  (  If6  (t):  t  ■  ; -t^.,  0;  J  is  an 
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open  cover  of  0)  and  8inc“  t*1**'  °)  u  c0***01*  *  “  finite 

autocover  [  M-  ...  ,  ^  U.)  )  of  l«V»  °1*  1  "  1#  " 

ti  tB 

let  Kt  •  {  ti  jt-tj  5  6t  i  t  t  0)  )•  Note  that  each  Kj  U 

«  coetpect  sutoeet  of  (-,0};  in  fact*  each  ^  ie  «  co«pect  aubeet  of 
[-twe,  0],  end  cleerly  u  Aic0'  lf  1  c  Ki  th<in 

*  i 

Hv(t)  -  $(^>8  <  ^5  •  tor  1  -  I.  ...,»,  let  UA  -  !  y:  y  €  E°» 

Hy  .  $<♦  )|j  <  Ji-  ),  and  let  A(K1#  denote  the  set  of  all  $  €  C 
*  Uk 

mapping  K.  into  U..  Then,  n  A(K£,  Uj  is  the  required  set: 
l  i  l»l 

Clearly  9  €  X  U^i  and,  It  1  €.  A  A^,  ^ )  then 

i*l  lmL 

||^(t)  -  $(t)||  <  2  l*^#)  f°r  t  in  1‘t^e,  0],  and  hence  p($,f)  ■ 
‘l*k 

£  m.  <  £  n*  +  2  onpl  IfcU)  "  -tk+l  -  t  *  'fcK  *  * 

k-cT*  k«k«  k-o 

^  ^  ♦  k#  aup(  ||p(t)  -  *(t)||i  -tk*  S‘^0)  < 

£  »,  ♦k*.2.|-»i|  s  I  s  >*/2*“/2 

k*K# 

re 

p,  or  p(qp,f )  <  p.  Thuo,  9  €  3l$.u)»  end  hence  A(Kt,  l^) 


c  S(cp iP )  • 


Q.E.D. 


APWKDIX  C 


In  this  appendix  ve  shall  prove  that  the  functionals  V  defined 


in  Equation*  (5-1)  and  (>5)  are  continuous  on  for  any  M  >  0, 
Theorem.  For  any  H  >  0,  if  V  is  defined  on  CH  by  either  Equation 
(5-1)  or  Equation  then  V  is  uniformly  continuous  on  Cjj. 

Proof:  Let  n  >  0  be  given.  Let  •  ry  «•  +  ...  n0,  Choose 


the  positive  integer  so  large  that 
2H2  t 

'“t1*  ■*  t€“  “a  k2“i  ' 


Choose  the  positive  integer  Ng  so  large  that 


2 

2H  T . 


bi  uTcjl)  lpelK+1  <  1  c  0  40(1  k^N2- 

Let  N  *  max  (N^  Ng) .  Let  T  -  max  {t^  :  tea)  and  let 


(*-!p,|  ) 


r— •  ^  j  \  I P 

n  '  2_  (i- jp, 

_  » 


(a&  ♦  bt|pj| )  +  bjij 


and  choose  6  >  0  such  that  5,  i  €  C„,  p(5,  5)  <  6 

n 

=>  sup (  115(c)  -  »(o):,  :  -  (2N+1)T  <  0  <  01  <  min 

This  is  possible  since  convergence  in  the  compact  open  topology  on 
(-«*,  0]  is  equivalent  tc  uniform  convergence  on  every  compact  subs*' 
of  ( -»,  0).  For  ever/  T,  i  t  CH  with  o(5,  i)  <  5  , 
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[V(i)  -  V(i)] 


1  rt. 


I  *‘(0)  p£<0  5  -  -  f^o)  tf(o) 


“  -Pkx  5  0 

-  }_J‘i  Z!'pt^k  /  *  bj  r<yt(a)cp(aj 

4«o'  k*o  0=-2{kl-l)T  t  *•  a«-t . 


dc 


5«0  k-o  "o=-2{ki-l)T  £ 

2(k+l}T  ^  , 

£  lpj|k+1  J  ^(oj^ojda 

'■  ~  cr=-(2(k+l)+l)t 


k^<3 


4 


£jttg  £^4^  f  *t(o)5(o)do  +  b J  f  ^o)i  o)do 
lefi  k*>  o— 2(  k+l)Tfi  a=-t 

-  - ( 2k+l )t ,  . 

^|p5|k+1  />(*)*(  a)do  j| 


k=o 


i  irt 


<  r  |^(0)  P$(0)  -  Vi 0)  P  $(0)[  + 


/  ®  -PkT , 

5«a\  k^0  a=-2(k+l)x 


U 

*tfa)^(o)|  da  +  b^[cpt(  a)^(  a)  -  vt( a)l( a)]  da 


00  -(2k+l  )x  ( 

b4  £  lp£  |k+1  f  [$*(<*)$( a)  -  d1' o)$(a)jcia 


k=o 


<  I  l9t(°)  P  <jp(0)  -  $l(0)  P  v(0)  +  9l(0)  P  $(0)  -  ^(0)  P  f(0)| 

-2kT. 


a  =-(2(k+l)+lJt 


4 


£  [a4  £  lp^k  J  I  IM a ) I!"  -  li^ o) |j2  | (la 

£e8  k=o  y  f  ir4-i  w 


a=-^{  k+J  ;t  : 


+  /|  l!<p(o)|!2  -  ||f(oH!t  !-a 


no 


<*  ||<f(0),  r($(o)  -  i(0)  I  ♦  ||<  pt(^(o)  -  1(0)  ),  *(o)>  I 

N-l  -2Xt^ 

+  E  f*{  E  Ipjf*  /  *  ^(o)l|do 

t«a  '  k=o  tj=*2(k+i)t^ 

O  N-l  -(2*+1)t5 

+  b  f  2H-(l9(fl)  -  $(o)||da  +  b  V  jp,  |k+1  f  2H-||9(c)  -  f(a)|jdtf 
a=-Tt  cr=-(2(k+l)+X)xt 

«  -2kT^ 

+  r  I  “{  e  k^/1  ^(°)|12  • ||f(a)|i?  id° 

tea  L  k=N  yo=-2(k+l)Tj 

*  -(2k+l)t, 

+  bt  E  |pt<k+1  / 1  ^(o)i|2  -  ^(o)ii2  ld0 

k=N  a=-(2(k+l)+l)tt  1 

<1  1^(0 ) || • -IIPII -mo)  -  |(0)||  +  |  1^11-119(0)  -  *(0)||-||*(0)|| 


1We  use  here  the  In-.. qua.  ivy  |  jj^?(  a )  ||2  -  |j$(a)j!2_l  <  2H.|j9(c)  -  f(c)||, 
vhich  follovs  immediately  frcm  the  inequalities  I  ll<p(  oTII  -  ||«(o)||  I  < 

(|cp(, c )  -  f(0)||  and  j  (|<p(o)!!  *  l|v(c)||  j  <  2H- 
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s  i  (2  f)  * 


3*,U  -  ip,jN) 

“ir~CT) "  ui +  b 


,I»||)  *  Vi] 


*  Zl4  +  rf~! a  t 

t«a  8  8 


£  £  f 

r  *•  r  +  r  =  e 

4  4  2 


That  is*  9,  f  €  Ch/  p(fj  , ;  <  f. ]v($)  -  V( v)  j  <  e. 

Thue,  V  la  uniformly  continuous  on  C„-  Q.E.D. 


Theorem.  For  any  H  >  0  If  V  in  defined  on  C„  by  either  Equation  (5-1) 

n 

or  Equation  (5-5)>  the?-  there  exists  a  real  number  t  such  that 
|  V  ( <p)  |  <  t  for  every  if  r  C5J. 

/  i  r-  +  b*!Pt  I)  ]  \  2 

Hoof:  Let  i  .  |||P||  ♦  2.  ^  '  lb  [f  '  +  Yl  * 


Then ,  if  q>  €  C„, 


«  -2KT| 

iv(«p)f  <2  HpH *  IMpC°  > II2  4  V  U|  ^  1 P4  Ik  f  ||^(  o)li2cio 


*  -o 


0— 2(K+1)t. 


o  fl0  -(2k+l)tj 

J Hv>( a) ||2dcr  +  t.  ^  1  of  ! V '  j  llqKtjJifdc 
o  -r^  *  -o  o=-(2(K+i  )+l  )t  t 


<  r  !'P| 


1)2  *  iK  i *  v/  *  k;  e  ^i'^S"2 


tea  k=o 


( ? ,Pil  *  E  [ (rrtfn  v.  *  tt^0)  ]) 


t-  3-E.D. 
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Electrical  networks  consisting  of  lumped  linear  and  memoryless  non-linear  elements 
and  an  arbitrary  number  of  lossless  transmission  lines  are  considered.  It  is  shown 
that  a  large  class  of  such  networks  may  be  described  by  a  system  of  functional 
differential  equations  having  the  form 

7ft)  -  T  (7  ) 

where  the  state  of  the  system  at  time  t^O  is  represented  by  7t,  a  noint  in  the 
snace  Cj|  (  f-**t03  ,  En  )  of  bounded  continuous  functions  manning  the  interval 

rrj  into  E0  t  with  the  comnact  open  topology,  and  the  function  f  manning  Cjj 

(  -••,  Oj  ,  E°  )  into  E°  is  continuous  and  locally  Linschitizian.  A  Lyapunov 
functional  is  presented  and  used  to  obtain  several  theorems  concerning  the  stability 
and  instability  of  the  equilibrium  solution,  7  ■  AT,  of  the  network.  Several 
exanrlcs  of  the  theory  are  presented. 
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